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Strong Turbulence by Renormalization 
I. KINETIC THEORY OF TURBULENCE- 

The hydrodynamical equations of turbulent motions are inhomogeneous 
and nonlinear in their inertia and force terms, and will generate a 
hierarchy. Dr, C. M. Tchen has developed a kinetic method to transform 
the hydrodynamic equations into a master equation governing the velocity 
distribution f(t,x,v) , as a function of the time t, the position x, and 

>- w ^ 

the velocity v as an independent variable. The master equation presents 
the advantage of being homogeneous and having fewer nonlinear terms, and 
is therefore simpler for the investigation of closure. 

After the closure by means of a cascade-scaling procedure, the 
kinetic equation is derived and possesses a memory which represents the 
non-Markovian character of turbulence. The said kinetic equation is 
transformed back to the hydrodynamical form to yield an energy balance 
in the cascade form, i.e. with a mode-transfer that is explicit. Note 
that the underlying mechanism is a "collisionless damping", similar to 
the "nonlinear Landau damping" in plasma turbulence. 

The theory analyzes the eddy transports along a gradient (i,e. normal 
transports), counter the gradient or without a gradient (i,e. anomalous 
transports) . 

The kinetic theory is described irt PAPER #1 for the incompressible 
turbulence, in PAPER #3 for the compressible turbulence, and in PAPER #2 
and #5 for the plasma turbulence. 

The normal and anomalous transports are described in PAPER #4. 
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II. APPLICATIONS 

Theory of turbulence developed above by Dr. Tchen served as a 
statistical foundation to the following problems. The application to 
sound generation is given in PAPER 6. 

The application to the propagation of light in a non-frozen turbulence 
is given in PAPER #7. In this connection two questions are of particular 
importance. The first deals with the remote sensing of the cross-wind 
by using two laser beams. The theory of remote sensing in non-frozen turbule 
derives a formula that can measure both the wind velocity and the eddy 
viscosity, while the formula based upon the hypothesis of frozen 
turbulence (Taylor's hypothesis), as well known in the literature, can 
only measure the wind velocity, but not the characteristics of turbulence. 
The second question deals with the frequency spectrum of phase fluctuations. 
The theory based upon the assumption of frozen turbulence (i.e. with an 
unperturbed path of light) gives a spectrum of -8/3 power, as well known 
in the literature. Our new result as obtained by treating the closure 
without the hypothesis of frozen turbulence (i.e. with a perturbed 
path of light) yields a spectrum of -2 power law, which is in a better 
agreement with available experimental data. 

Other applications to an atmospheric medium with particulates are 


also given in PAPER #7. 
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MECANIQUE DES FLUIDES. — Modele ciitetiqiie dc turbulence dans un jluidc 
incompressible. Note (*) de Chan M. Tchen, transmise par M. Marcel Barrere. 


Nous developpons un traitemcnt slatistique pour I’etude de la turbulence dans un fluide incompressible 
obeissant aux equations de Navier-Stokes, cn considcrant la pression comme un potentiel d’interaclion entre 
les elements de fluide. Un proccde d’echelles divise une fluctuation en trois rangs representant les trois 
processus de transport : evolution macroscopique, propriete de transport et relaxation. Nous obtenons une 
fermeture par relaxation, et deduisons une equation cinetique pour la fluctuation de rang macroscopique de 
la function de distribution. La solution donne la fonction de transfert et la viscosite tourbillonnaire. Comme 
application dans la zone d’inertie du spectre d’energie, nous trouvons la loi de Kolmogoroff et son coefficient 
numerique. 

]Ve develop a statistical approach for iiwestii/atinfi incompressible Navicr-Stokes ttirhiilcnce by consiilcriiifi 
pressure as a potential of interaction between fluid elements. A scaling procedure divides a fluctuation into 
three ranks, representing three transport proce.sses: macro.scopic evolution, transport property and 
relaxation. tVe obtain a closure by relaxation, and derive a kinetic eguation for the distribution function of 
macroscopic rank. The solution yields a transfer function and an eddy viscosity. Application to inertia 
subrange finds the Kolmogoroff law with numerical coefficient. 


I. Introduction. — La representation cinetique de la turbulence dans un fluide 
incompressible est basee sur I’hypothese suivant laquelle la pression agit comme un potentiel 
d’interaction entre deux elements de fluide, par analogie avec le potentiel entre deux 
particules dans la theorie cinetique d’un systeme discret. Cette representation conduit a une 
hierarchie de fonctions de distribution a plusieurs points de I’espace, semblable a la hierarchie 
de BBGK (Bogoliubov, Born, Green, Kirkwood) a plusieurs particules [('), (‘‘)]. 

Dans la theorie dc Bogoliubov, la fermeture de la hierarchic consistc a degenerer une 
distribution d’ordre superieur en distributions d’ordres inferieurs dont la justification peut 
etre donnee par un dcvcloppement en puissances d’un petit parametre. Lc succes de la theorie 
cinetique de la turbulence dans un systeme discret est, en grande partie, du a I'existence de ce 
parametre qui n’existe pas dans notre systeme conlinu. Pour cette raison, nous developpons 
une theorie dc fluctuations, par un procede d’echelle des fluctuations de fonctions de 
distribution dont nous chcrchons la fermeture par un processus de relaxation asymptotique. 

Comme la fonction de distribution normalisee a I’unite est equivalente a la densite de 
probabilite, ce modele cinetique coincide avec le modele probabiliste. 

II. Representation ciNfeTiQUE. — Les equations de Navier-Stokes et de continuite, dans 
un fluide incompressible avec vitesse u, pression p, densite p et viscosite cinematique v, 
peuvent etre considerees comme dcs moments de I’equation microscopique : 

(1) (d, -l-L)N(t, X, v) = 0 avec L = v-V-1-E -t)„ — vV^ 


La fonction de distribution N(t, x, v)s 5 [v — u(t, x)] satisfait a la condition de 
normalisation 


( 2 ) 


i/vN (1, X, v) = 1 . 

J - no 


Comme lc champ est defini par 

(3) E(f, x)=p-‘ Vp(t, x)=-ry(x/x', v'){N(t, x', v'}, 
on peut ecrire (1) sous la forme 

(4) (5,-t-v -V-vV^jNft, X, v) = d„g(x/x', v'){ N(t, x, v)N(t, x', v')}. 
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ou 


g(x/x', v'){. . .} = 


4n 


ti xVv'V|x-x'|-‘ V'V' : v'v'{. . .} 


est un operateur integral, et d, = d/d,, V = d/d\, V' = d/d\', d„=d/ov. 

Dans un mouvement turbulent, la distribution N =/ +/ se decompose en une distribution 
moyenne/et une fluctuation /, ou/ peut s’appeler densite de probabilite, au titre de la 
normalisation a I’unite d’apres (2). 

En prenant la moyenne de (4), au moyen de I’operateur A , et cn deduisant la fluctuation au 
moyen de I’operateur A = l — A, nous obtenons : 

(5 a) (d,+AL)7=-Ad„-(E/), 

(5 b) (5,+AL)7=-E-5„7. 


On dispose alors de deux methodes : la premiere mcthode consiste a trailer la distribution 
moyenne en ecrivant (5 a) sous la forme 

(6) (d, + AL)7=a„g(x/x', v'){</(t, x,.v)/(r, x', v')>}, 

qui engendre une hierarchic. Cette forme a ete aussi dcduite par Monin (*) et Lundgrcn (^) 
par d’autres moycns; la seconde methodc que nous envisageons ici consiste a trailer les 
fluctuations par un procede d’echelles, en ecrivant : 

(7) /=/“>' + /', /' = /<*>+/”. 

Les trois rangs/*®’,/' et/" representent les trois processus de transport en plusieurs 
cascades : evolution macroscopique, propriete de transport et relaxation, ayant leurs durees 
de correlation 

( 8 ) 

d’apres la discrimination par les operateurs A‘°’, A', A". Les echelles instantanees des rangs 
peuvent se depasser, mais leurs echelles statistiques (8) sont distinctes, de telle fa9on que les 
intensites 


<E(0)2> 


= 2| dk'S(k'), 


<u“’>^> = 2 


dk'F(k') 


sont reparties dans un domaine (0, k) des spectres S(/c') et F(k'). 

Par I’intermediaire des operateurs A‘°', A', A" appliques a (5 b). nous obtenons les 
equations des rangs, et en particulier, par relimination du rang microscopique caracterise 
par A', nous obtenons I’equation cinetique suivante : 


(9 a) (5,+A<°'L)/‘'»=-E<o>-3„7+A<°>d„-0'.{d„/‘°'}, 


avec 


(9 b) &' = 


dr<E'(r, x)U(r, r-T)E'(t-x)>. 
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Cette equation cinetique gouverne la fluctuation et differe de I’equation cinetique 
suivante, deduite de la meme maniere 

_ _ ao 

(10) (5, + AL)/=d„.0.{d„/} avec ®'= dt< E(t, x) U(f, t-i)E(t-T) >. 

J 0 

Cette equation peut etre consideree comme une forme fermee de la hierarchic (6). Une 
equation cinetique pour la distribution / 2 (t, x, v, x', v') a deux points peut aussi se deduire 
d’une maniere analogue. Elle pourrait servir a developper une fonction de correlation et par 
consequent, une fonction spectrale. L’equation cinetique de en (9 a) a I’avantage de 
se dispenser de et de permettre, directement, la determination de la fonction spectrale. 

III. Turbulence isotrope et homog£ne. - Une correlation lagrangienne 
<E'(t, x)U(t, t—i) E'(r — x)> apparait en (9 b). Son calcul suit une trajectoire qui est 
perturbee par la turbulence deE(f, x). Dans 1 = 1 +1‘‘” + 1“’ + 1", nous distinguons le 
parcours au rang submicroscopique l"(t) qui traduit une diffusion des trajectoires, fournit 
une relaxation ainsi que les autres parcours qui correspondent a la convection. Par 
I’isolement des rangs E' et 1 "( t ), I’expression (9 b) peut etre decouplee et s’ecrire sous la 
forme scalaire qui est la trace du tenseur isotrope 

*1 -* 00 CC'^ 

(11) D'=2 (lx dco'dk'S(o)'. k') 

Jo jj 

X <exp-ik-l"(T)><exp [-vk'^ T + i(co'-k' • v)x-ik' ■(T + l‘°' + l<“)] >. 

la moyenne se calcule par une probabilite de transition obeissant a une equation de Fokker- 
Planck, et S(co', k') est la fonction spectrale dans I’cspace de (co'. k'). Une turbulence faible 
permet de negligcr la relaxation, tandis qu’unc turbulence forte permet de negliger la 
convection. Dans le dernier cas, nous reduisons (1 1) a 


(12) D' = |r(^^j ”dk'S(/c')/coM^'). 


* 1 "l‘/3 


La resolution de I’equation cinetique (9 a) et la methode des moments donnent alors le 
bilan d’energie 

(13) = 

La fonction de transfert se compose d’un produit de viscosite tourbillonnaire 

(14) K'=|r(^^jr0) *dk'S(/c')/co[,^ 
et d’une fonction de vorticite 

(15) R<°> = 2 rdk'/c'2F(fc'). 


Si nous remplapons (13) par le bilan d’energie qui egale au taux de dissipation e dans la 
zone d’inertie, nous trouvons la fonction spectrale 

(16) F = 
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j conforme a la loi de Kolmogoroff (*). Le coefficient numerique A pent etre calcule par une 

equation qui relie E a u en partant de (3) et cgalement S(/c) a F (k). Le calcul donne ainsi la 
, valeur du coefficient A = 1,6, pour k/ k^= 100, ou fc*’ est I’cchelle interne de Kolmogoroff 

I gouvernant le diametre des plus petits tourbillons avant leur dissipation par viscosite. 

Ce travail est fait sous 1’ auspice de National Aeronautics and Space Administration. 

I (*) Seance du It septembre 1978. 

(‘) A. S. Monin, Appl. Math. Mech., 31, 1967, p. 1057. 

(^) T. S, Lundgren, Phys. Fluids, 10, 1967, p. 969. 

I C ) F. R. UuNlCH ct B. Ya. Lyubimov, Zh. Eksp. Tear. Fiz., 55, 1968, p. 951; Sov. Phys.-J .E.T.P., 28, 1969, 

) p. 492. 

(*) V. M. Ievlev, Dokl Akad. Nauk S.S.S.R., 208, 1973, p. 1044; Sov. Phys. Dokl., 18, 1973, p. 117. 

-A (’) A. N. Kolmogoroff, C. R. Acad. Sc. U.R.SS.. 30, 1941, p. 301. 

J Graduate Canter and City College, 

City University of Neu> York, New York, N. Y. 10031, U.S.A. 
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PHYSIQUE MATHEMATIQUE. — Fermeture par relaxation d'un modele cinetique de 
plasma en turbulence forte. Note (*) de Chan Mou Tchen, presentee par M. Andre 
Lichnerowicz. 



La decomposition en cascade repetee permet de deriver une equation cinetique de Fokker- Planck pour un 
plasma turbulent. Le calcul du propagateur et la solution de I’equation cinetique determinent les coefficients 
de transports (diffusivite et viscosite tourbillonnaires), moyennant une fermeture basee sur un processus de 
relaxation responsable de I’approche 4 I’equilibre. Ensuite la loi spectrale k~^ est pr6vue, et gouverne le 
couplage des fluctuations de vitesse et de champ electrostatique. 

A scaling procedure by repeated-cascade derives a kinetic equation in the form of Fokker-Planck equation for 
turbulent plasmas. The calculation of the propagator and the solution of the kinetic equation determine 
transport coefficients(eddy diffusivity and eddy viscosity), with a closure based upon a relaxation process that is 
responsible for the approach to equilibrium. A spectral law k~^ is found to govern the coupling between 
velocity and electrostatic field fluctuations, and is in agreement with experiments and numerical modeling. The 
same kinetic method may serve in determining the k~^ spectrum for velocity fluctuation in compressible 
turbulence and the Kolmogoroff law k"*'’ in incompressible turbulence. 

I. Cascade r6p6t6e. — La fonction de distribution / {t, x, v) de vitesse v dans un plasma 
sans collision est gouvernee par I’equation de Vlasov : 

(1) (d,+v.V + E-d)/=0, Vsd/dx, d = d/dy. 

Le champ self-consistant E qui sera precise plus tard, est suppose fluctuer sans valeur 
moyenne. 

On decompose la fonction de distribution en 

(2) /(r, X, v)=/(t, v)+7(t, X. v), 

ou / est la moyenne et / est une fluctuation, et le champ fluctuant en un macro-champ E‘°’, un 
micro-champ E', et un champ submicroscopique E", comme suit : 


r 

a 


f. 



(3) E = E‘“' + E', E' = E<‘> + E", (E<®»-1-E<‘> = E,). 

La discrimination de ces composantes ou rangs se fait par les moyennes d’ensemble, 
representees par les operateurs fibres A“”. A', A", qui s’applique aussi bien a/(t, x, v). 
L’application des operateurs filtres A*®’, A' a (1) conduit au systeme suivant ; 

(4) (a,-t-A“»L)/<°>=-E“»-a/-5.<E'/'>. (d,-FA'L)f'=-E'-d{f-hf°>), 

ou/= +/', et I’elimination de/' nous donne I’equation cinetique sous la forme de Fokker- 
Planck ; 

(5) {d, -1- A<°> L) /‘°> = - E«» • d /-H A‘°> d • D' • a /“>' ,L = vV-f-E-d, 
avec une diffusitivite 

(6) D'= dt < E'{l, x) U(t, t-T) E'(f-T)>, U = propagateur. 

J 0 

La decomposition en cascade repetee permet de definir une fonction spectrale S{k) 
d’energie electrostatique, telle que 

'k 

(7) <E'°'"> = 2 dk'S(k'), 

J 0 
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appartient aux nombrcs d’ondc plus pclits que k. landis quc la diffusivile D' est effectuee par 
le micro-champ E' aux nombres d’onde plus grands quc k. Finalemcnt une relaxation 
controle le temps necessaire aux proprieles de transport pour s’approcher a I’equilibre; elle 
regit les nombres d’onde plus grands que k' et par consequence sera gouvernee par tin champ 
submicroscopique E". Les nombres d’onde fcet l:'sont priscomme variables independantes. 
Les rangs E“”. E', E" sont autorises a outrepasser Icurs echcllcs instantanees. mais leurs 
durecs de correlation doivent ctre nettement separees comme echelles 

statistiques. Une fonction spectrale F (k) pour la fluctuation de vitesse u pent etre defmie 
d’une manicre analogue. Le coefficient % en (7) provient de la transformation de Fourier 
bornee. 

11 faut remarquer que les decompositions (2) et (3) sont usuelles en turbulence 
hydrodynamique : (2) s’appelle la decomposition de Reynolds, et (3) s’appelle la 
decomposition en cascade repetee qui est une generalisation de la cascade de Heisenberg (*). 
Les trois rangs caracterisent les trois processus dans I’ordre d’echclies dccroissantes ou dans 
I’ordre croissant au titre alcatoire : 1’evolution macroscopique, le developpement des 
proprietes de transport, et la relaxation, forment un couplage qui rejoint la hierarchic 
cinctique de B.B.G.K. (^). 

II. Fermeture. — D’apres la theorie de diffusion, la diffusivite devient asymptotique, 
i.e. D'(r>x') = D'(r -» oo), quand f>x', cc qui est vjai puisque t est rechelle d’evolution 
de en (5). Par le meme raisonnement une diffusivite qui possede un temps disponible x 
plus petit doit choisir le rang D"(x) pour devenir asymptotique, i.e. 

D"(x>x;')=D"(x-oo), 

correspondant au rang l"(x) du parcours l(x) = l, (x)-f l"(x). 

Le calcul du propagateur U transforme la diffusivite (7) sous la forme de Fourier : 

(8) D'= ' dx ” dm'dk'x <E'(co', k') E'(-co', -k')> 

• 0 » * “ "TO 

X < exp-i k’- l"(x)>exp i[o)'x-k' -(vx-t-l,)], 

qui contient une convection libre exp i[o)'x — k'*(vx-t-I,)] et une relaxation 
< exp — / k'-l" (x) > qui pent se decoupler du champ E' a cause de leur difference de rangs. 
L'effct dc convection libre se manifesto en turbulence faible, comme il est d’usage dans la 
theorie quasilineaire, et I’effet de relaxation domine en turbulence forte, ce que nous traitons 
id. La moyenne de la fonction exponentielle est calculce par une probabilite de transition 
gouvernee par une autre equation de Fokker-Planck semblable a (5). Le calcul 
transforme (8) en une equation integrale dont la solution donne la trace du tenseur isotrope 
suivante : 


( 9 ) 



dk' S(k') 



3/4 


111. Structure SPECTRALE. — Considcrons un plasma quasi-neutre. compose d'ions d(^ 
masse M, au scin d’clectrons chauds dc temperature T_ ayant une densite ii qui est rclicc au 
champ E= — Vvl/f suivant la loi d’cquilibre n = ;io exp(vl//Co), oii Co=(>cT_/M)'^^ est la 
vitesse thermique. 

On peut defmir une fortetion d’enstrophic pour le potcntiel par I'exprcssion 
J“’' = M < (V vl/'“’)^ >. qui sc reduit a J quand k s’accroit indefiniment, et une fonction 
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d'enstrophie analogue R'°’ pour la vitesse. Lcs equations d’evolution d’energie cinetique ct 
potentiellc peuvent se deduire comme moments de I’equation (5) de Fokker-PIanck 
gouvernant la distribution d’oii I’on ctablit les equations du bilan d’energies. En 
omettant les details, on obtient dans le cas de stationarite statistique et dans la zone 
d’echange : 

(10) R<°> = Cte, C'®>-PrK' J‘“' = Cte, 

oii la viscosite tourbillonnaire est ; 

et la fonction d’echange est : = M D*°’. 

L’energie cinetique est excitee par les fluctuations clcctrostatiques occasionnant un 
cchange O'” entre Ic spectre d’energie cinetique ct Ic spectre d’energie potentielle. Les deux 
energies poursuivent leurs dissipations tourbillonnaircs par Icurs transferts K' R"” 
et Pr K' J’®’ dus au termc d'intcraction non lincaire dcs ondes a travers leurs spectres 
individuels. Le transfert dans le spectre d’energic potentiellc s’effcctue avec une propriete de 
transport modifiee par un nombre de Prandtl turbulent Pr. 

II apparait done que Ic spectre d’energie cinetique sc trouve dans le regime de 
dcvcloppemcnt par Ic gain d’echangc oil R*®’ s 0, tandis que Ic spectre d’energie potentiellc 
tend vers son regime terminal par lcs deux pertes qu’ellc subit, ou J'®’ s J a grand k. Le flux 
d’energie, decrit par la forme differentielle de (10), peut s’ecrire, par approximation ; 

(12) R'®* = 0, -C<®^-Pr K' J = 0, 

oil (■)= (l/(lk. L’enstrophic J dcfinit un puits dans Ic processus d’echange cl de dissipations 
tourbillonnaircs des deux spectres, et remplacc le taux de dissipation dans la iheorie de 
turbulence hydrodynamique (®). Dcs equations (12), on tire les solutions 

(13) F(k) = PrJk-\ S(k) = Cte (J/iIq M)^ 

La loi /c“® est bien rcconnue dans les mesures (“*) ct calculs numeriques (^). Dans la 
comparaison avec lcs calculs numeriques, il faut remarquer que la loi ~ obtenue 
numeriquenient comme composante de scrie de Fourier conformc a la prediction (13). 

Notre methode cinetique peut servir a determiner la loi spectrale et la loi de 
Kolmogoroff dans les cas de turbulence compressible et incompressible, 

respectivement. 

Ce travail est fait sous I’Auspice de National Aeronautics and Space Administration. 


(*) Seance du 13 fevrier 1978. 

(‘) W. Heisi-.nuero, Z. Phyaik, 124, 1948, p. 628. 

(^) J. Salmon, Ann. Inst. Henri Poincare, 27. 1977, p. 73. 

(^) A. N. Kolmogoroi r. C. R. AcaJ. Sc. U.R.SS.. 30. 1941, p. 301. 

(*) C. M. Toie.n, H. L. PixsLLi ct S. E. Larsen, Ris0 Report, n" 365, Research Establishment Risp, DK-4000 
Roskilde, Denmark. 

(’) J. J. Thomson, R. J. Fap.iil, W. L. Kruer ct S. Bodner. P/iys. FtniJs. 17, 1974, p. 973. 
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A KINETIC APPROACH TO COMPRESSIBLE TURBULENCE 

C. M. Tchen 

City College, City University of New York, N.Y. 10031, USA 

We develop a kinetic theory of strong turbulence, by 
decomposing fluctuation of distribution function into 
scaled components, or ranks, to represent three transport 
processes: macroscopic evolution, transport property and 
relaxation. Closure is obtained by formulating a relaxation 
which determines the approach of eddy diffusivity to 
equilibrium. We solve the kinetic equation and derive tran- 
sport functions of turbulence in adiabatic gas: two trans- 
fer functions describe turbulent mode transfer, and a 
coupling function describes the build-up of kinetic energy 
at the expense of internal energy. Their balance in the 
inertial subrange of homogeneous and isotropic turbulence 
in an adiabatic gas yields a k”^ law for kinetic energy and 
internal energy. Their intensities predict a k“^ law, in 
agreement with numerical computations from one dimensional 
model of compressible turbulence. 

I. INTRODUCTION 

By considering pressure as responsible for the interac- 
tion between fluid elements. Monin , Lundgren2, Ulinich, 
Lyubimov^, and Ievlev^ have described a kinetic treatment 
of turbulence. The kinetic hierarchy generated is similar 
to the BBGK hierarchy in plasmas, and has to be closed by 
degenerating a high order distribution function into lower 
order ones. 

The success in this direction will depend on finding 
a small parameter, as in the case of plasmas. However, 
this parameter is not present in neutral turbulence. In 
order to avoid this difficulty, we develop a fluctuation 
theory, by scaling a fluctuation into components of various 
degree of random character, or ranks, to represent macro- 
scopic evolution of distribution fluctuation, eddy diffusi- 
vity from microscopic fluctuation, and eddy relaxation from 
submicroscopic fluctuation. The relaxation takes the form 
of two streamings (free and dynamical) and submicroscopic 
diffusion which is made asymptotic for the approach of eddy 
diffusivity to equilibrium. Such an asymptotic diffusion 
renders the relaxation frequency deterministic and degene- 
rates a high order correlation into lower order ones, thus 
achieving a closure. 

This closure produces a kinetic equation with an expli- 
cit diffusivity. By solving the kinetic equation, we can 
calculate the transport functions consisting of a coupling 
function and two transfer functions, governing the energy 
balance for homogeneous and isotropic turbulence in an 
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adiabatic gas. Eddy viscosity is found. We derive spectral 
distributions of kinetic energy and internal energy, from 
velocity and density fluctuations, respectively. The k"^ 
law of these energies agree with a one-dimensional numeri- 
cal modeling developed by Tokunaga.6 

II. FLUCTUATION THEORY OF TRANSPORT 

A. Interaction Between Fluid Elements 

We introduce a distribution function 

f(t,x,y) = f(t,x) S[v-u(t,x)J (1) 

of velocity v at time t and position x. It gives density 
and momentum by *“ 

Jdv f = f , /dv y f = |>u . (2) 

A self-consistent field 

E = VP = - •^^//d^’dv’^C|x-x’( ,x* ,v*)f (t,x’y») 

— p .. ^ ^ I «, w ^ — ~(3) 

from pressure gradient Vp can be written. Here the operator 

[x-x’( ,x’ ,v*) = I x-x’l (3^+y’ V* *V’ (4) 

represents the interaction between two fluid elements in 
compressible fluid. It will reduce to | x-xM "^(v’ 
in incompressible fluid, and to the Coulomb potential 
|x-x’| in a plasma. It is easy to verify that the Liouvi- 
lle equation 

(V+L)f = 0, with -L = vV+E-"^ (5) 

will reproduce the equations of continuity and momentiam by 
means of moments. Our transport theory, which will derive 
eddy transport functions and coefficients, can neglect 
viscous effects. Here "i. ='b/dt, V = "^/^x, . 

T- w w ^ ^ 

B . Scaling 
We can write 

E = E + E , ^ = E + E » , E ’ = E E" . (6) 

Vv* 

The components E ^ , E’ and E" partake in macroscopic 

evolution, eddy diffusivity, and relaxation, respectively. 

The various components , or ranks , may overlap in 
instantaneous scales, but have statistically distinct 
correlation times, or scales, in the decreasing order 

%' > T” . ( 7 ) 

signifying their increasing degree of random character, 
and, consequently, the quasi-stationarity of one macrosco- 
pic rank with respect to the next random rank. 
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The ensemble averages A ^ , A increa- 

sing finer scales help in discerning the_rai^s E(°^ E’, 
E”, by means of filter operators 2 = 1-A, A^°^= A -A, 

A’ = 1-A , A"= 1-A^ . . ° 

o’ 1 

C. Kinetic Equation 

By the application of A^°\ A^ , we transform Eq. (5) 

into 

t ^ j 

= - E».;b(f+fW) +A^<Lf^> , 

and by the elimination of f * , we derive a kinetic ec[uation 

+A^°^L)f f + A^°^D’l^^f (9) 

L w ^ jtr-ir 

with a diffusivty 

D» = j^^dT A^<Je’ C t)A»U(t,t-T)E» (t-T))> (10) 

which may also serve as an integral operator on functions 
within brackets ■[•••}. The propagator U indicates a 
Lagrangian representation. 

III. CLOSURE BY RELAXATION 
We can write Eq. (10) in Fourier form: 

D* =j^'*^dT^ d4J*dk’^ A^{e’ (- o»,-k’)EUo’,k*)Ajh^, 

n i") 

where h is a kernel, and ^ is a factor of truncation ^ ■' 

in the Fourier transformation of a function which is 
quasi-stationary within a certain time interval or space 
interval, or both. 

We have the following kernels: 

A.H = A,hh , A,h = A,h"h, , 

1 1 V’ 1 1 1’ „ -ik».^» 

h = exp i(tb’-k^*v) , h = expT-k’ (T)?,A^h"=A^ e ^ 

( 12 ) 

The total kernel h consists of an Eulerian kernel h 
representing a streaming of velocity v, and of a Lagrangi- 
an kernel h from the perturbed pathJ^=X+^". In its turn, 
the Lagrangian kernel h is separated* into“h, and A,h”, 
representing a dynamical streaming of path and a 
relaxation by diffusion of the random path '"Jt”, 
respectively. We shall rely upon this micro-diffusion for 
obtaining the approach of D* to equilibrium and closure. 

From the inequalities 



T(," < X < <r t , (13) 

the upper limit t can be replace by oo in Eq. (11) . For 
the same reason, while x is larger enough for E" and A^^h" 
to reach asymptotic diffusion, i.e. od , it is 
allowing Ej^ and h^^ to afford a dynamical streaming only. 

Since ^ 

kjh = h^, for weak turbulence (14a) 

- A^h" , for strong turbulence , (14b) 

an interpolation formula 

Ajh = Aj^h"h^ (14c) 

can be written by neglecting dynamic streaming. 


The relaxation A^^h" can be calculated, by means of a 
Fokker-Planck equation for the transition probability 
P(t,£*^), with a diffusion consistent to the same process 
as governed by Eq. (11) . We find: 

Aj^h" = exp(-«jj”^T^) , with | k’^D". (15) 

The relaxation frequency and the asymptotic diffusivi- 

ty 

D" = tr y'“dx»^du”dky A^<(E"(-^l’^-k’0E;'(<J’^k" A 2 h’" 

= tr dr y dk" ^ A^(e" (-k") E» (k"))> exp (- ^ T ^) (16) 

are found independent of v in strong turbulence. 

By introducing a spectral function F^(k’), defined 
by the trace 

t^dk";rA^(E"(-k")E”(k» )) = 1/ dk" Fj^(k") , (17) 

and by performing the x’ “ii^legration, we reduce Eq. (16) 
to 


= r 8 r/4 \ f 
9' {3J J 

|_ 3 ■Vt 


dk" Fg(k")/(|k"^)^^^ 


]3/. 


(18) 


Note that 


1 E" I ^ in strong turbulence 
/V/ I E" I , in weak turbulence . 


( 19 ) 
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IV. PHYSICAL BASIS OF BALANCE BETWEEN KINETIC ENERGY 
AND INTERNAL ENERGY 

Introduce an auxiliary variable 

« ““h “o = fot’ (20) 

having the dimension of velocity, in an aidabatic gas, with 


an equation of state 

P = Pjf/f/ 

= ^ - ( 21 ) 

The equation of continuity 

can be written as 

(3^+u*7)w = - w V-u , (23) 

to form the equation of energy 
h \ +V^(u^%fw^j= -V-(up) u-E . (24) 

We can verify that 

%w^ = c^T (25a) 

is the internal energy , and- that 

%Yw^ = CpT = ^ (25b) 

is the enthalpy playing the role of an effective potential, 
such that 1 

E = - ^ Vp = - Vf . (26) 

Hence Eq. (24) will be called equation of internal energy. 

The two forms taken by Eqs. (23) and (24) are analo- 
gous to the equations of momentum 

(3^+u-V)u = E (27) 

and of kinetic energy 

+ 7j(u^-Vu^) = U'E . (28) 

In homogeneous turbulence, the two energy equations 
reduce to 

•it (29a) 

■it Su = C , (29b) 
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where 

C=(/>u-E^ (29c) 

is the coupling function which describes the build-up of 
kinetic energy '8 = % u^"^ at the expense of internal 


energy . Note that 




(29d) 


The energy balance described above is written in the 
non-scaled from, i.e. k = od . The corresponding scaled 
form, i.e. k = finite, must retain internal transfers 


\ 2 
s, 


as f o3JLows : 


^ - T 

u 

(30a) 

= - C - T . 

w 

(30b) 


(o) 


(o) 


The transfer functions T , _ 
mode-transfers. We have rhe conditions: 


govern the internal 


C^°^= C, 


TuC°)= 0, 


T 0, at k=cD , (31) 

w 


so that the system of Eqs. (30) reverts to Eqs. (29) at 
k = 00 , as expected. 

We shall write the kinetic equation (9) in the form: 


(*i^+L)f^°^ = = Jj 

where ^ 

5 - E^°^.J)f 
SO that the solution 


(o) 


f^°^t,x,v) = d i;’y^dcj^dk* J 


+ (32) 

D»{^ ^f (33) 


(o) 



{cd't-k'-ii) 


= fg^°^ + fp^°^ (31^) 

consists of two parts, as contributed by J ^°^and J , 
respectively. ^ 

The presence of h and of the T -integration facilitate 
the transformation of - ±cJ into 

m(x’) = 2 S' (-[:*) - ik’*v (35) 

from time differentiation, obtaining 
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C'tjX,v)=A^°y Jp ,k* ,v) 


X m(x,’)Aj^he 


-i(<y^t-k* *x) 


C36) 


In terms of distribution functions, the transport 
functions in Eqs. (29) - (31) can now be written as 

- hz ''i* 

- /dv k 

(,(o) = A /* dv V . (37c) 


V. DERIVATION OF TRANSPORT FUNCTIONS AND TRANSPORT 
COEFFICIENTS 

The procedure of calculations of transport functions 
requires a closure. This is made by the use of the inter- 
polation formula 

A h = A h»h (38) 

O O V ^ 


on an analogous basis as Eq. (14c). The partial integra- 
tions with respect to v and the integration with respect 
to time are made. The details of calculations will be 
omitted. 


For strong turbulence, we find the following transport 
functions : 

Tu^°^ =y (39a) 

T = Pr“^ -> 7 ’ R (39b) 

w / w 

C^°^ = 3 f D^°^ , (39c) 


the transport coefficients: 

7’ = I [r Cj)] F^(k*)/ 6 )j 3 "^ (40a) 


(40b) 



- 8 - 


the vorticity functions : 


R zf dk’ k’^F (k’) 

u ^ 

R zf dk» k»^F (k*), R = R 

w Jg w ^ ’ w w [k=oo ’ 

the spectral functions: 

^ dk» F^Ck’) = A dv 

/•k=CD P 

/g dk* F^Ck’) = % 

and the turbulent Prandtl number 
Pr = 16/3 Y - 


(^la) 

(i+Oh) 


(42a) 

(42b) 

(43) 


VI. ENERGY BALANCE AND SPECTRAL STRUCTURE 

In statistical equilibrivun, the inertia subrange in 
compressible turbulence is governed by the following 
coupled system 


q(o) _ T (o) = c _ £ 

u 

-C - T = -C , 

w . ’ 


(44a) 

(44b) 


for kinetic energy and internal energy. The system is a 
generalization of the energy balance 


T £ 

u 


(o) 


in incompressible turbulence, for which C and C vanish, 
where f is the rate of energy dissipation. 

The system describes a build-up of kinetic energy at the 
expense of internal energy for, an amount C^ , to be 
followed by mode-transfers T ^ "^and T in respective 
spectra. In differential formv we have: 

C^°^ - ^’R = 0 (45a) 

- y u / u 

-C^°^ - Pr“^^'R^*^°^ - Pr"^^*R^^°^ = 0. (45b) 

By this coupling, the tv\?o spectra in the inertia subrange 
are in two different regimes: the F^ spectrum is in the 
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regime of production where its vorticity function R is 
negligible , while the F spectrum, after having suffered 
two losses, falls in the regime of eddy dissipation where 
its vorticity function becomes asymptotic: R ^ R . Under 
these circumstances, the energy balance reduces to: 


= 0 

- Pr"^i5»R = 0. 

/ w 


(46a) 

(46b) 


The upper dot denotes a differentiation with respect to k. 

With the aid of Eqs. (39) and (40) , we can solve the 
system of Eqs. (46). The solutions can be expressed in 
terms of Pr"'^R or in terms of the constant relaxation 
frequency ^ _ ? r „ u. ^ 

as a parameter. 

We find the spectral functions 

-2 


= =>uf 

r- _ -3 

^E ^E ^ ’ 


5 a 

’ u 


'= 2-52 


(48) 


ap = 9/ r (V3) = 10.08, 


the intensities 


Eu’ = ^ =>ur 


<1'^'^ =>£ 


.■t ,,-2 


(49) 


and the transport coefficients 

o,’ = a tJ * k“^ 
y y D > 

D» = 3 k"^ . 


a = 4r(4/3) =3.57 

7 


(50) 


We can also calculate the intensity of enthalpy 
fluctuation 


from Eqs. 


(26) and (48) , and the internal energy 



2 -2 
EJ ’ k ^ 
D 



v-1 
3.17 / 


(51) 


( 52 ) 
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from Eq. (25b) . 

_2 

The k law for kinetic and internal energies in Eqs. 

(M-9) and (52) are in agreement with results of numerical 
computations for one-dimensional model by Tokunaga.° 

VII. CONCLUSIONS 

Instead of closing a kinetic hierarchy of mean distri- 
butions, we develop a fluctuation theory, and close a 
hierarchy of correlations. The closure is facilitated by a 
scaling procedure into ranks. When the microscopic diffusion 
of path becomes asymptotic, a relaxation is obtained, pro- 
viding equilibrium to diffusivity and rendering it deter- 
ministic, and in so doing, it closes the hierarchy. The 
kinetic equation of macroscopic distribution is then solved 
for deriving transport functions and transport coefficients. 
The balance between kinetic energy and internal energy 
derives the spectra. In a customary kinetic treatment 
without scaling, one would require a doublet distribution 
function to calculate a velocity correlation, and a Fourier 
transformation to derive the spectrum. The present method 
by scaling fluctuations, a singlet distribution function 
suffices for the same purpose. 

With the transport coefficients determined in this 
manner, the method could also lend to analyse the distri- 
bution of velocities around a mean value and the deviation 
from normality. This problem will be left for a later 
opportunity. 
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TURBULENCE. - Transports turbulents normaux et anormaux dans un milieu stratifie. 
Note (*) de Tchen Chan-Mou, presentee par Alexandre Favre. 

L’etude ost consacree a la derivation analytique des flux de transport dans un milieu stratifie. La memoire, le 
couplage entre Ics fluctuations de vitesse et de temperature et les proprietes tensoriclles apportent aux transports 
anormaux un caractcre contre-gradient ou sans gradient, outre Ic transport normal du type Boussinesq suivant un 
gradient. 

This Note presents an analytical derivation of the fluxes of transport in a stratified medium. The memory, the 
coupling between velocity and temperature fluctuations arid the tensorial properties of transport give rise to an 
anomalous counter-gradient or without gradient transport, in addition to the normal transport of the Boussinesq type 
along a gradient. 

I. Introduction. — Les transports anormaux de quantile de mouvement ([1] a [4]) et de 
chaleur qui nc suivent pas la relation de Boussinesq onl etc mis en evidence dans des 
ecoulcmcnts turbulents cn laboratoire [5], dans ralmosphcre ([6], [7]) cl cn astrophysique [8]. 
On pent expliquer ce phenomene par un effet de memoire, un couplage entre les fluctuations 
de vitesse et de temperature, et les proprietes tensoriclles du transport, que nous examinons 
separement dans la suite. 

II. Memoire ET coui'LAGii. — On distingue les valeurs moyennes u, T, de vitesse cl de 
temperature, et Iciir fluctuations turbuicntes u, 0. Introduisons tine vitesse tiuxiliairc 
uiscG/Tq, un champ E=— V(p/p), deux gradients y = Vn) + (r//c)2, L = ViI) + 2.oZ. 
et deux operateurs differentiels et L„ tels que los A(u + u). V w, L„usL„,u — E, oil 
cs(RTq)*^^ est la vitesse de propagation, R cst la constante des gaz parfaits. To’ 
correspond au coefficient d’expansion thermique, tj est I’acccleralion de la pesanteur; 
2.0 est un taux adiabatique, z = (0, 0, 1) cst Ic vcctcur unite, p cst la masse volumique 
moyenne, p cst la fluctuation de pression, et A cst I’operatcur de fluctuation. Les equations 
fondamentales de quantile de mouvement et de chaleur peuvent s’ecrire sous la forme 
suivante ; 

(1) (i), + L„)u= -u.Vu + ii'Yr ((', + -u.L. 

Par integration de (1), nous obtenons les vitessesu, ic, au moyen desquelles nous calculons 
les flux _ 

( 2 ) 

<Uj«;>= -^(K ;#(. 


Les coefficients de transports 


(3) 


1 ^ u, u> 

(M„); = 


dx<K^(t, x)U„. „(f, f-T)H,(f-T)>; 


dT<n;(f,'x)U„(l, t-T)u!(t-i)> 


agissenl comme operateurs. L’inlegration sc fait suivant les Irajccloires qui soul peiTurbes 
par les champs turbulents, comme I’indiquent les propagaleurs U„, qui sonl les inverses de 
L„, Ainsi la forme lagrangienne, qui porle sur les gradients cn (2), represente un effet de 
memoire. Cel effet a servi de base pour expliquer le dccalagc [9] entre flux et gradient. Le 
couplage entre vitesse et temperature el les proprietes tensoriclles donneront d’autres 
transports anormaux. 
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III. Proprietes TENSORiELLES DE TRANSPORT. — Lcs cocfficiciits clc transport (3) sont 
exprimes avec ou pour un temps fmi qui est Ic temps d’approche a 

I’equilibre, appele temps de relaxation ou o),;;*, en introduisant lcs notations 

Qj,(co, k) = x<"j(“- -k)>, RjSx<i/j((o, k)ti'{-to. -k)>, 


oil le facteur x provient de la transformation de Fourier bornee. Nous pouvons exprimer (3) 
sous la forme de Fourier [10] : 


(4) (K„.J,,= 


dcodkco“i„Q;,(cp, k). 


(M„),.= 


CO 

deodkeo" ‘ Rj(co, k), 

- CO 


et ainsi transformer (2) en 

Qji(co, k)= k)V,Ui + (o;‘ Rj(to, k)Yi, 

5 

Rj((i), k)= Qy,(o). k)3.,. 


(5) 


Ce systeme peut sc resoudre en fonction dc la composantc diagonale Qjj pour donner 
(6) Qji=-Qjj^jir. = aj.v(3.) 


et par consequent pour obtenir : 

(7) (K„. Jj,.= -(K„. {M„),.= cj, ,>(?.)}. 

avec Nous remarquons que lcs coefficients dc transport (K„ „,)jj operent sur 


} 

i 

] 


j 


(8 a) 

Sj,r=( r„ r,, - r„)/( r., r„ - r,. r,> 

et 


(8 b). 

I (X) S kj-Sji, Xi - S J>i Xr. 

I a Ji, ( V M j ) = ( V J— S j-,v V i - S jVi V J u , , 

avec 


(8 c) 

= 5si + Vj Ui/0)„ + X,y ,/co„ CO 


Une theorie d’echelle pour la fermeture par relaxation [10] permet de determiner la 
structure des frequences de relaxation to„, co„. Ici il suffit de considerer (8 c) comme un 
gradient non dimensionnel donne, sans entrer dans la structure analytique des frequences de 
relaxation. 

IV. Transports normaux et anormaux. — Par la substitution de (7) en (2), nous 
aboutissons aux flux 

Rappelonsque(K„ sont les coefficients de transport sous forme dc Fourier (4) et operent 
sur lcs gradients compo.scs (8 b), pour donner un transport normal dii type Boussinesq avec 
memoirc et les transports anormaux. Lc transport de quantitc de mouvement est plus 
complique que celui de la chaleur, et peut devenir plus facilement anormal. Si de tcls 
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transports anormaux predominent, les Ilux peuvent changer dc signe, de sorte que la 
viscosite effective pent formellement devenir negative [8]. 

Nous illustrons ccs effets par un excmple, ou la presence des gradients non nuls V3U1, 
V| It’, et X,3 = V3tt’+^o reduisent (9) au systemc suivant ; 


( 10 ) 


<«3«1 > = -(KJ33{A, '(V3U,+tJ„X,3)}; 

<W3U;> = -(K,„)33{ Ar'(>.3-r/„ V3U,)}, 


aveca„_ Vj u)/co„ Aj = 1 + auCi,t,. Les flux verticaux(10)acontre-gradient resultent d’un 
gradient horizontal Vi w qui apparait dans a„;„. 

Dans une stratification uniquement verticale avec gradients V3 u , et V3 lu, les transports 
deviennent 

(11) ^ “1 )= ~(l^u)33 { 1^3 tt 1 } ; ^ U 3 U) > = — (K „,)33 { V 3 U'+ A-o } , 

appartenant au type de Boussinesq avec memoire. L’absence de tous gradients V u , V lu, sauf 
le taux adiabatique, peut encore donner lieu a un flux de chaleur 

(12) <“3'i’> = -(K»)33 ^o- 

Si les gradients varient tres lentement, le systeme (1 1) se reduit aux transports normaux sans 
memoire 

(13) <«3»1>=-(KJ33V3«„ <«3«’> = -(KJ33(V3>I>+-^o). 

qui sert a donner la justification de la theoric de similitude de Monin et Obukhoff [11]. 
Cependant cette theorie de similitude ne traduit pas la dcpendance quantitative des 
coefficients de transport a la longueur de stabilite. Notre theorie d'echclie [10] permet de 
determiner la structure des coefficients de transport et de la relaxation pour analyser cette 
dependance. 

Ce travail a etc fail sous les auspices dc la National Aeronautics and Space Administration. 
Resume d’un texte qui sera conserve 5 ans par les ArchTves de I’Academic cl donl copie peut etre obtenue. 


(•) Remise le 25 juin 1979, acccptec le 17 dccembre 1979. 
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Abstract — An investigation of the spectral structure of turbulence in a plasma confined by a strong 
homogeneous magnetic field was made by means of a fluid description. The turbulent spectrum is 
divided into subranges. Mean gradients of velocity and density excite turbulent motions, and govern 
the production subrange. The spectra of velocity and potential fluctuations interact in the coupling 
subrange, and the energy is transferred along the spectrum in the inertia subrange. Applying the 
method of cascade decomposition, the spectral laws k~^, k~^, k~^ are obtained for the velocity 
fluctuations, and k~^, k~^, k~^^^ for the potential fluctuations in the production, coupling and inertia 
subranges, respectively. The coefficient of Bohm diffusion is reproduced, and its role in electrostatic 
coupling is derived. Comparison is made with measured power laws reported in the literature, from 
Q-devices, hot-cathode reflex arc, Stellarator, Zeta discharge, ionospheric plasmas, and auroral plasma 
turbulence. 

1. INTRODUCTION 

Collective anomalous properties of turbulent plasmas are found important in 
characterising plasma motions in fusion devices, in astrophysics and in ionospheric 
heating. Observations have been reported on the spectral structure of turbulence 
for density, electric potential and electrostatic field fluctuations, (e.g. Bol, 1964; 
Chen 1965; D’Angelo and Enriques, 1966; Berkl and Grieger, 1967; Reece 
Roth 1971; Robinson and Rusbridge, 1971; Leonard and Linnerud, 1972; 
Gurnett and Frank, 1977). It is our purpose to present a hydrodynamic 
formulation of the theory of turbulence for comparison with these observations. A 



kinetic formulation was given previously by Tchen (1978a, b). 

Theories of weak plasma turbulence have been frequently reported in the 
literature (Kadomtsev, 1965). Theories of strong turbulence (Dupree, 1966; 
Weinstock, 1969) are limited to weak couplings. For strong couplings, a proce- 
dure was developed for smoothing the fluctuations by a repeated-cascade as an 
advanced Reynolds’ decomposition (Tchen, 1969). This procedure was applied to 
the Navier-Stokes model for a two-dimensional magnetized plasma (Tchen, 
1976). Here we shall generalize the problem to three dimensions and deal with 
the spectral structure of the fluctuating components in the plane perpendicular to 
the magnetic field. Further details of the analysis presented in the following 
sections are available in the form of a report (Tchen et al., 1977). 

2. MODEL EQUATIONS AND CASCADE DECOMPOSITIONS 

The motion of ions, of density n, velocity v, temperature Tj, mass M, charge e, 
un4«r an external magnetic field B = (0, 0, B) and a self-consistent field 

E=-V4i, (1) 

* Temporary address; Dept, of Phys. and Astronomy, University of Iowa, Iowa City, Iowa 52242, 
U.S.A. 
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is governed by the following equations of conservation of mass and momentum: 


— + V • (nv) = 0, 
dt 


nM^^+v • Vv^ — en(— V<|) + vxB)= — KTjVn. 


(2a) 

(2b) 


The hot electrons of temperature T, can maintain quasi-neutrality by flowing 
along the magnetic line of force, according to the equation 

(3) 


(4) 


(5) 


(6a) 
(6b) 

The system of conservation equations (2) and its reduced form (6) will serve as 
a basis for formulating a turbulent transport theory, as being separate from 
collisions or molecular dissipations. Therefore these dissipative effects are not 
written in (2) and (6), but will, however, be restituted as sinks in energy balance 
later on. 

If the plasma motion is strictly two-dimensional in the plane perpendicular to 
the magnetic field, we can reduce the system (6) into an equation for V x v of the 
Navier-Stokes type without pressure. Such an eqimtion can serve as a simple 
model of turbulence. However, if the above condition is not fulfilled, the complete 
system (6) has to be considered. This latter problem will be treated here. 

We decompose the velocity v into a mean velocity v and a turbulent fluctua- 
tion u, as 

v = v+u, (7a) 

and subsequently we redecompose the fluctuation into a macroscopic component 
u*®’ and a random fluctuation u', as: 

u = u“” + u'. (7b) 

The decomposition (7a) is known as Reynolds’ decomposition in turbulence, and 
the decomposition (7b) is an extension. The macroscopic component u'"’ and the 
random fluctuation u' can be screened by means of a distribution function, so that, 
after an ensemble average. 


kT, Vn = tie V(j) 

for a Boltzmann distribution. By introducing 

i/» = a In (n/no) 

with a normalizing constant iiq and 

a = ^K(T, + Ti)IM 

we can rewrite (2a) and (2b) in the form: 

Sip „ „ 

— + \‘Vip = — aV • V, 
dt 

Sv „ „ c „ 

— l-v • Vv = -aVip-\ vxB. 

dt M 
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|((«;)^>= [”dfc'F(k'), 

Jk 


represent two consecutive portions of a spectral function F(k') for velocity 
fluctuations separated by a wavenumber k considered as an independent variable 
ranging from 0 to Before this ensemble average, and u' may have 
overlapping length scales, but should predominantly lie in separate spectral 
portions. The same decomposition applies to t|r, so that 

t|/ = til + t/t<°) + (//'. (9) 


m^°Y)-Cdk'G{k'), dk'G(k') (10) 

Jo Jk 

are two consecutive portions of the spectral function G(k') for i/r-fluctuations. 
Similarity, for E-fluctuations, we have: 

Ke/M)WV>= f'dk'GE(k'), Ue/My((E;f)= f dk'G^(k’). (11) 

Jo Jfc 

Decomposing the system (6) into macroscopic and random components, the 
following equations are obtained in the form of a cascade for components in the 
plane perpendicular to the magnetic field: 

Diul”'- (12a) 

- adju<°^-(u/djip'} (12b) 

(a. + Vidj)Ui' - (o^Siiu/ = ^ Ei' - u/d,{Vi + uf”) (13a) 

(d, + Vjdj)ip' = -adiU/-u/d,(if/ + tl/^“^) (13b) 


The fluxes 


D, = a, + (v + u<°>)-V, a,=- 

01 




-E = -aViJ/. 
M 


<n/n/>, (u/^'> 


are Reynolds’ stresses accounting for the statistical effects of random fluctuations 
upon the evolution of macroscopic components and 1 /^*°’ in (12). In turn, the 
random fluctuations are produced by the gradients of macroscopic variables, as in 




I 
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(13). The discrimination between the two systems (12) and (13) is obtained by a 
scaled ensemble average. As the fluxes (15) are of macroscopic rank, they are 
omitted in (13), governing fluctuations of random rank. 

3. TRANSPORT COEFFICIENTS 

The system (13) can be integrated formally to give the random components: 

Mi' = -[5s(Ui + Uj“’)] dTt/„(f, t-T)u/(t-T) +■£:[ dTl/„(f, t-T)£j'(t-T) (16a) 

Jo Jvz Jq 

t^' = -[3j(i/?+i/'"”)] dTi;(t, t-T)u/(t-T)-a3i[ .dTU(t, /-t)U|'(i-t) (16b) 

Jo Jo 


which serve to form the gradients of Reynolds’ stresses: 

(u/a,.u/> = -K,,'i3;a,(0i + O (17a) 

'a, a,(.i;+ ,/,<“>) (17b) 

and the cross-correlation 

-<£/n/)-D,,'. (17c) 

The eddy viscosity in x-space, the diffusivity in x-space, and the diffusivity in 
volocity space, which are assumed homogeneous, are: 

dT{u;(t,x)UJt,t-r)u/U-r)) (18a) 

A) 

\j,'= dT<u/(f, x)l/(f, t-T)u/(l-T)) ( 18 b) 

Jo 

D,/ = (^)'|^'dT<£/(l,x)t;,,(U-T)£/(l-T)), ( 18 c) 


respectively. The fluxes (u/u/) and <u/i/x') in (17a) and (17b) are calculated 
approximately by only retaining the gradient of the macroscopic quantities 
transported, i.e. a,(ui + uf’) and aj(t// + 1 //‘“’), while the cross correlation (17c) is 
calculated approximately as a transport of non-gradient type. The propagator 
t-r) implies a Lagrangian representation of the correlation function follow- 
ing a trajectory that is perturbed by u in a magnetic field. For a strong magnetic 
field, we can write 

f~T)~exp[-(wfc +iwjT], (19) 

where = eB/M is the cyclotron frequency, and is a relaxation frequency 
governing the approach to equilibrium of the transport property, and it is 
approximately 

- «c- (20) 

Thus we can write (18a) and (18c) in the form of traces as 


I K'=ak(Wj'(/,x)]^)w,, ' 

1 and / e 

D' = aD(^) 

1 


( 21 ) 

(22a) 



Strong turbulence in low-(3 plasmas 


821 



respectively. Using the definition (14), we can rewrite (22a) as 




(22b) 

where 





(23) 


is found as the Bohm diffusion. The transport coefficients K^°\ and D“”, 
which are defined from macroscopic fluctuations and E^°\ can be written 
similarly. The numerical coefficients and are not determined. It should be 
remarked that most analytical theories of turbulence have not reached the stage 
where numerical coefficients can be determined. The present analytical descrip- 
tion of plasma turbulence, as different from similarity theories which are based 
upon dimensional arguments, aims at a description of turbulent transport proces- 
ses and the derivation of the spectral structure of the transport coefficients. 


4. ENERGY BALANCE 

Multiplying (12a) and (12b) by and and averaging, we find the 
following equations of energy balance in collisionless form: 



(24a) 


(24b) 

These equations consist of the production functions 


1 = 

(25a) 

J P»'> = A<“T/, 

(25b) 

a coupling function and the transfer functions 


rj” = 

(26a) 

= A'J<"*. 

(26b) 

j The gradient and vorticity functions arc: 



(27a) 

] R‘'” = (a,MlT) = 2[ dk'fc'^F(k') 

(27b) 

= = dk'k'^Gik') 

Jo 

(27c) 

respectively, with asymptotic values: 


^ k 

R=((a.u.)Z) = 2 dk'fc'"F(k') 

(27d) 

j^(0.,/,)2> = 2['' dk'k’^G(k'). 
Jo 

(27e) 


The transfer functions (26), which arises from nonlinear mode-couplings, 
represent an energy cascade, as was postulated by Kolmogoroff (1941) and 
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Heisenberg (1948). Transfer functions in the form (26) has a kinetic basis (Tchen, 
1978a, b). They do not exist in theories of weak turbulence. 

All the transport functions are considered homogeneous in the locally 
homogeneous plasma. 

5. SPECTRAL STRUCTURE 

A. Classification of spectral subranges 

We shall distinguish between universal and non-universal ranges in a spectrum 
and shall be exclusively concerned with the universal range. The non-universal 
range depends on the conditions for the particular experiments. 

The universal range can be subdivided into production, coupling, and inertia 
subranges. The individual subranges are investigated separately. 

B. The production subrange 

Mean gradients of velocity and potential will feed energy into the fluctuations 
for further transfer across the individual spectra. In order to describe this 
transport explicitly, we write the energy equations in the following differential 
form, where the upper dot represents differentiation with respect to k; 

= 0 (28a) 

= 0 (28b) 

obtained by retaining the production and transfer functions in (24). 

As the production occurs at low wavenumbers, we may neglect the last terms 
in (28a) and (28b) on account of small and J*°’. The simplified equations then 
become 

rjK^°' = K’R‘°\ (29a) 


yielding the solutions 




F(k) = const 
G(K) = const 


C. The coupling subrange 

Density fluctuations which appear as potential energy can excite the kinetic 
energy and give rise to a transfer down the spectrum. This circumstance defines 
the coupling subrange. The governing transport functions in this subrange are thus 
the coupling function and the transfer function. From (24), we have, in differential 
form 

D‘o)-t).°> = 0 (32a) 

-D"”-T«’> = 0 (32b) 

or, with the substitution of (22), (25) and (26), 

A,j‘"*-JC'R‘“'-K'R‘°> = 0, (33a) 

= (33b) 
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Because the i/» spectrum is dissipated at a rate controlled by the Bohm diffusion, 
the following approximations can be introduced 

^(0)^0, and A'«A„ 


reducing (33) to 

A7<°>-K'R<°> = 0 

(34a) 


-A7"’>-A'J = 0 

(34b) 

or, after addition. 

-K'F<°>-A'J = 0, 

(35) 

Substitution of (21) into (35) gives 



A' 

(36) 

with 

(Oq^ = 

(37) 


As the left-hand side of equation (36) is a function of F and k alone, we 
obtain the solution 

F(k) = const (38) 


On the other hand, equation (34a) gives, by means of the definition (21), 

2k"‘G(k) = const — F<“>. (39) 

AcWc 

Substituting the solution (38) for F(k), we find from equation (39): 

G(k) = const (<uo“/A^oj,.)k"* (40) 

corresponding to a spectrum for the electrostatic field fluctuations 

Gg(k) = const (a^ajo“/Ac<u<.)k"^ (41) 

using relation (14). ~ 

We can conclude that the t/r-spectrum (40) falls off much faster than the 
M-spectrum (38) in the present subrange, on account of the effective dissipation by 
the Bohm diffusion in the t/r-spcctrum alone. The velocity and the electric field 
spectra share a k~^ law, as evidenced by available data. 

D. Inertia subrange 

The inertia subrange is characterized by a constant transfer of energy across 
the spectrum, i.e. 

T^«’=eu (42a) 

T*;' = (42b) 



X'F“”=£o 

A7‘"> = £.,. 



or, with the definitions (26) 


(43a) 

(43b) 
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Here e„ and are the rates of energy dissipation, or sinks. They serve as 
parameters governing inertia and dissipation subranges. 

The rate of energy dissipation can be identified as 

e^^D = KJ (44) 

from (22b) and (27c). Note that is not due to coliisional losses but to a 
collisionless Bohm diffusion. This is to be expected, since the equation of 
continuity (2a) is conserved. On the other hand, the equation for the momentum 
(2b) may have a collision term giving rise to e„. The particular form of the 
collision term (kinematic viscosity by self-collisions, or damping by collisions 
between species) is not of interest here. 

The energy balance derived in (43) has a form similar to that postulated in the 
theories of Kolmogoroff (1941) and Heisenberg (1948), using dimensional 
considerations. Our transport theory, which derives from (43) and the transfer 
functions, is based upon a kinetic foundation (Tchen, 1978a, b). Since the 
structure of our eddy transport coefficients differ from that assumed by Heisen- 
berg (1948), we expect different solutions of the equations of energy balance 
(43). 

To obtain the solutions, we write (21) in spectral form 


and rewrite (43a) as 


Jk 

f dk'F(k')2{ dk"k"^F(k") = e„ 

Jk Jo 

|-k i 

dk'k'^F(k') = const CuO),.-— ; . 

“ dk'F(k') 

■ k 


After differentiation with respect to k, we obtain 
k^F(k) = const ~ 


• V’/ p fOO ^2 

dk'F(fc') 

■ h J 

• ao 

dk'F(k') = const (e„w,.)''^k~‘. 

J|c 


Differentiation of (48) finally yields 

F(k) = const (e„w,.)''^k"^. (49) 

On substituting (49) into (43b), we obtain 

G(k) = const (50) 

As the ijj spectrum is quickly dissipated by Bohm diffusion, as mentioned in 
Subsection 5C, we may expect that the inertia subrange (50) is not fully developed 
under certain circumstances. 
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6. COMPARISON WITH EXPERIMENTS 

We shall compare our theoretical predictions with data on turbulent spectra 
reported in the literature. In most experiments, the spectra are measured as a 
function of frequency rather than wavenumber. This introduces the problem of 
Eulerian-Lagrangian transformation. If the streaming velocity predominates, we 
obtain a linear relation between the spectral functions in w and k spaces, 
according to Taylor’s hypothesis, or frozen turbulence, which is very widely 
accepted in the interpretation of turbulence measurements in plasmas and fluids 
(e.g. Kofoed-Hansen and Wandel, 1967). The streaming velocity consists effec- 
tively of the mean velocity v and the velocity representing that portion of the 
spectrum of larger scale than the subrange concerned. Consequently, the Taylor 
hypothesis is valid for plasmas with a strong drift v, or, in its absence, for 
subranges (coupling and inertia subranges) of large wavenumbers. The hypothesis 
may not be valid for the production subrange (of small k) without a sufficient drift 
V. Under circumstances where the Taylor hypothesis is not valid, a nonlinear 
dispersion will prevail and this gives the frequency spectrum as a Fourier 
transform of the Lagrangian correlation in (18c); 

Gn(o)) = ;^(^)'|^ dre-‘-HE,'(t,x)UJt,t-r)E/(t-T)}. (51) 

For w = 0, this spectrum becomes linearly proportional to diffusivity: 

GE(tu = 0) = -D“”(k = oo). . (52) 

TT 

We note that Ge(w = 0)>0, while G^(k = 0) = 0. This feature is observed in all 
experimental data. 

An analytical study of the nonlinear dispersion relation requires the calcula- 
tion of the propagator U„ and of the Lagrangian correlation (51). Although the 
propagator has been investigated (Tchen, 1978a, b), the calculation of the Lan- 
grangian correlation has not yet been made. However, a comparison between the 
theoretical prediction of the spectrum in k-space and the available data in 
w-space can be justifiably made applying the Taylor hypothesis for the coupling 
subrange, with spectra 

F(k)~k-^ GE(k)~k-^ G(k)~k-^ (53) 

according to (38), (41) and (40), respectively. The predictions are found in 
reasonably good accord with experiments. 

In a strong streaming which is the necessary condition for the Taylor 
hypothesis, it is not clear what correction the non-uniformity of the stream may 
bring to the conversion between co and k representations. It may conjectured that 
a weak non-uniformity may not be important in the conversion of small scale 
turbulence. 

It is to be noted that, in the lack of a nonlinear dispersion relation, the Taylor 
hypothesis of frozen turbulence is widely accepted, not only in comparison 
between theory and experiment relating to w and k-spaces, but also in turbulent 
propagation of electromagnetic waves, in remote optical sensing of wind velocity, 
and in laser doppler measurement of fluid velocity. 
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Fig. 1. — Turbulent spectra obtained in Q-devices. All curves show the amplitude of the 
fluctuations as a function of frequency. The curves reproduced from D’Angelo and 
Enriques (1971) are marked Vand O: those of Berkl and Grieger (1966) by +. The 
measurements marked V and O refer to fluctuations in floating potential, while the 
experiment denoted + measured fluctuations in the ion saturation current to a Langmuir 
probe. The spectral index for the amplitude variations is denoted p, i.e. A(f)~ The 
amplitude is obtained by taking the rms output from a constant bandwidth filter. 



Fig. 2. — Turbulent power spectra obtained in a hot-cathode reflex arc (( ) by Chen 

(1965) and in the Etude stellarator ( ) by Bol (1964). The measurements refer to 

density fluctuations. The spectral index for the power spectrum is denoted u, i.e., 

,P(/) 





; 
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The inertia subrange may not be fully developed in experiments for the reason 
mentioned in Subsection 5C, and therefore the spectra (49) and (50) will not be 
compared with experiments. On the other hand, strong mean gradients arc 
unlikely to be found in plasma devices, except e.g. in the ionosphere. Therefore 
the spectral predictions (30) and (31) in the production subrange may not be 
easily found in data from plasma devices. 

With the reservations stipulated above we compare our theoretical predictions 
with experiments. Figure 1 shows spectra obtained in conventional Q-devices 
reproduced from D’Angelo and Enriques (1966) Berkl and Grieger (1967). 
The reported spectra show the amplitude variation as a function of frequency. We 
note that the measurements of Berkl and Grieger (1967) have a large ExB 
rotation of the entire plasma column, rendering Taylor’s hypothesis effective. 

Figure 2 shows turbulent spectra for the density fluctuations obtained in a 
hot-cathode reflex arc (Chen, 1965) and the Etude stellarator (Bol, 1964). 
Spectral measurements of the fluctuating electric field in the Zeta discharge 
(Robinson and Rusbridge, 1971) are shown in Fig. 3. 

Figure 4 shows a spectral analysis of turbulent density fluctuations in barium 
plasmas released in the upper atmosphere (Leonard and Linnerud, 1972). 
Because the data were reduced by photographic technique, it was possible to 
present the spectrum as a function of wavenumber in this case. 

Finally, Fig. 5 shows satellite data for the fluctuating electric field on high- 
latitude auroral field lines (Gurnett and Frank, 1977). 

In view of the difficulties associated with frozen turbulence mentioned above, 
the lack of supporting experimental data for determining parameters, and because 
diverse conditions of turbulence can produce a k~^ law (e.g. in neutral compressi- 
ble turbulence, see Tchen (1979)) we must caution any overemphasis of the 
agreement between theoretical predictions and experimental spectra. 





10 ’ 


10 ^ 10 ^ 
f (kHz) 



Fig. 3. — Power spectrum for the fluctuating electric field in the Zeta discharge (Robin- 
son and Rushhridge, 1971). 
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Spectral Structure of Noise Generated by Turbulence* 

C. M. Tchen 

City College, City University of New York, New Yox'k, N.Y. 10031 

ABSTRACT 

From the Lighthill equation for the generation of soiuid by turbulence, 
we investigate the spectral structure of sound px'essure fluctuations 
in a turbulent velocity field tliat possesses Kolmogoi’off ’s spectrum 
F (k) k ^ The pressure spectruja is found to be f|Ck)/vk where 

k is wavenumber. For the titans format ion of k-sijectriim into c<)-spectrum, 
a nonlinear dispersion relation for non-frozen turbulence is derived. 

The structure of Jl(k) contains a contribution from large scale motions 
which are linearly dispersed by free streaming (i.e. frozen turbulence) 
and a contribution from small-scale motions which are nonllnearly 
dispersed by diffusion process. The theory of transformation yields 
cO , The result compax’es favox’ably v^/ith measurements. 


* 

This work is performed under a grant fx’om the National Aex'onautics and 
Space Administration NGS 1485. 



I. INTRODUCTION 


The theoretical foundation of turbulent jet noise has been laid 
in the v?ork of Lighthill h] . Concerning the turbulent structure of noise, 
a large body of experimental measurements of two point pressure 
correlation has been cax'ried out by Maestrello [^ J. An analysis of the 
correlation has been developed by Ribner [^3j . The theory used an 
empirical Gaussian correlation function for tlie turbulent velocity 
fluctuations and found a good agreement with experiments. Because of the 
difficulties of correlation measurements at small distances between tx'io 

points, spectral measurements in frequency have been performed . 

- 5/3 

From the Kolmogoroff law k in the turbulent velocity field, 

- 7/3 

dimensional reasonings can predict a spectrum k for tlie micro-fluctuatioi 
of sound pressure, where k is wavenumbex’. The dimensional theoi’y does not 
provide detailed dynamics of sound emission. 

Since theory deals in k-space and measurements are in O -space, a 
comparison requires a space-time transformation, l.e. a dispersion 
relation. A linear relation based upon frozen turbulence is well known^G ^ 
but is not valid for small scale turbulence. We slaall attempt to resolve 
these difficulties Jiere. 


II. INHOMOGENEOUS WAVE EQUATION 


In incompressible, liomogeneous and isoti.'opic tui'bulence , tlie 
pressure fluctuation p(t,x) is related to velocity fluctuation u(t,x) 


by the Poisson equation 


- V^p = f, ^ » 

(1) 

through a source 

r z W • 

(2) 



-3- 


Here P is the constant density in tlie turbulent mediiun. Equation (1) can 
be generalized to 

( 3 ) 

for describing the generation • of noise of infinitesimal amplitude by 
turbulence ^1 J . The linearized equation (3) has a constant speed of 
propagation c and a turbulent driving source r . Tlie problem of 
acoustic turbulence which is fully nonlineax’ is treated separately {.7 J . 
Beside these two extremes of linear sound emission, or self -noise, and of 
acoustic turbulence, there exist otliei* models, including sound scattering. 

A derivation of these models is given in APPENDIX, including a rederivation 
of (3) . 


The solution of the inliomogeneous wave equation (3) is 
or, in Fourier form. 


(4a) 


(4b) 


wliere _ 

Here and in the following, tlie limits of integration are understood to 
extend to the whole available domain, unless oL'horwa.se Indicated. 

The relation (4b) serves. as a fundamental basis for analysing the 
spectrum of sound pressure fluctuations, as emitted in an isotropic and 
homogeneous turbulent field, whose velocity fluctuation is known to possess 
the Kolmogoroff spectrum. 


III. SPECTML STRUCTURE IN FREQUENCY AND WAVENU^U3ER SPACES 
A. Source 

In tlris I .. M.,. . . T . 



-ll._ 


(23 in terms of the spectrum of the turbulent velocity field. To tliis end, 
we decompose the velocity u into a constant mean velocity u, a macroscopic 

(o) 

fluctuation U. and a microscopic fluctuation u’ , as 


(;0) / 

U. =r •#' 6t i~ Ll 

*•%,« 

The operators , A’ select the ranks , u* , and are 

A^°^ = 1 - r - A’ 


( 0 ) 


A' 


.=. 1 - A 


to be regulated by ensemble averages A<^*-*^, A^ of increasingly 

finer scales, as denoted by the ranks in the operators A, A^ . 

We can write the macroscopic source 

A‘'’U 

A/u'u'>] , 


the condition of incompressibility 

= 0 


V • u 


C7) 


(8) 


and the stress 


U 


H 


(0) 


/ <fuv> - - yv 

where ^ ^ is an eddy viscosity wliicli is deterministic and stationary in 
t,x, to sliow tliat 

)~0>^ 


and to i>educe (7) ,to 


(o3 (o3 (o3 

r ^ u ^ u ^ • 

In wavenumber space, this is 

- - f ^k'(krkp^Tc<!p%yU-^'J 

./u, / - } ^ 


(10a) 


(10 b) 



We can calculate 




Simplifications will be made by using the following properties: 

(i) isotropy 

(ii) homogeneity 

4(u.%i ic‘tA"j> =• X a(u.%') U.^°^t-k'J>^Ck'i-k‘'J . 

^ ^ ^ VV 

(iii) incomxjressibility of turbulence 

k • u.^°>(h ) = O 

(iv) quasi-normality 

,/,^J (®A/ c°j\ / K/, 


Here ■ 


7, -Wxf 


( 11 ) 


(12a) 


(12b) 


(12c) 


(12d) 


(13a) 


is a factor of truncation of Fourier transforma Irion of a function wliich 
is s tationar'y within an interval of space 2X. When bhe Fourier trunsformat i 


is extended also, to time within an intei'val 2T , the factor concerned is 


A f (j/ 


(13b) 


instead of (13a) . 
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We introduce an energy spectrum F (k) j such tliat 

J = folk' A<\{°Kk')f> 

= ^ f^dy F( k'J , (Wa) 

'‘o 

with 

F(k) = 2-7Tk^^ A <^)u^°^ (k) ) A (14b) 


in the inertia subrange of isotropic turbulence, and A = l.G. It ±■ollov^;s 
an enstrophy 'function- 

=|”dk» k’^^ A <fju^°^ (k’)| 

rk 2 

= 2 / dl<’ k' F(k») . (15) 

''O 

Similarly, we introduce a spectrum S (k) for the fluctuations of tlie 
source 

A = I'dky^ A (]<’) I 

= 2 / dk’ S(k*) , (Iba) 

''o 

and a pressure spectrum HOO » such that 

A = /dk’ A(|p^°\k»)) 

^k 

= 2 cU<»ri(k’) . (Ibb) 

U 

The quadruple correlation in (11) cun be reduced by the use of 
(12a) - (12d) . It 'is shown that tlie first product 

<[u^°^(k’)u^°^ (k-k’)) (-k")u^^„^, 

■ in (12d) gives no contribution, and the remaining products yield 





4 


In the notation of (15) - and (10a) , we can rewrite (17) in tlie form: 


_ 2 . 


K 


(18) 


giving, by differentiation. 


S(k) = - F(k) . 


This relation is equivalent to 




or 


(19) 


( 20 ) 


= f 'X, ^jl '> ' 

B. Pressure 

From (Mb) and (21) , we derive tlie speeLral structure of pressure 




r-/u (o 


= i f ?!:°’ k'"' ^Hkc) K A , 

^ 'o Z 


(22a) 


1 
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and its intensity 


=. ff'pM' i°>(kj //<!«' %'W/!>. 


(22b) 


IV. DISPERSION RELATION 
A. Eulerian and Lagrangian Representations 
We shall, investigate the dispersion relation, i.e. a relation 
between frequency and wavenumber in turbulence, by considering tlie 
following transform 


(23) 


betvv^een the Lagrangian correlation at two tines t-TT ^ind t on tlie left 
hand side and the Eulerian spectrum ^A^ju (^^l<*)j ^/through a 
kernel on the right hand side. The propagator U(t,t-T) gives a 
Lagrangian form to the velocity u^°^ (t-T) which follows a perturbed 
trajectory. In Fourier space, the propagator becomes a Lagrangian kernel 


'll = A, ]i’ li 
1 u 


(^'l.a) 


consisting of: (a) an Eulerian kernel 


, _ -ik’*u 
li = e ~ 
u > 


(2Hb) 


v;hich represents a free streaming (detei’ministic) of a frozen turbulent 

pattern by a mean velocity U , and (b) a relaxation kernel 

' 3 3 

A^h» = e D ^ , 


(24c) 



which represents a diffusive perturbation of the trajectory at a 
relaxation frequency 

dependent on the diffusivity D". In the inertia subrange of strong 
tui’bulence, tliis diffusivity is determined to be 

D" = c < 

‘-D t , 

where £ is the rate of energy dissipation, and Cj^ = ^.2. *5 


(25a) 


(25b) 


is a numeric an coefficient. In weak turbulence, the kernel (2Mc) should be 

i. 

replaced by another dynamical kernel. This dynamical kei’nel is negligible 
in the present strong turbulence. The calculations of tlie propagator and 
the kernels f m the basis of closure by relaxation. For specific 
details, see Ref. 8. 

B. Diffusion Time 

The Lagrangian correlation, (23) defines a frequency spectrum 




i/COT 




(20a) 


whicli becomes 


with the use of (2Ma) • Here the dispersion time is defined by the 


(20b) 


Fourier transform of h(x) , as 


■=: ^ fcl-U ’^(x) -e, 

= ^fc{z k) t] + 


ccoT 


1 
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t>y (24) , and satisfies the normalization condition 


h(cj,k) ^ I • 


( 28 ) 


It may not be too wrong to replace Ijy 

estimate (27) by means of the following approximate form: 


1 






J> 

This function may be considered as a distribution of k wliich peaks 

s 

at a value of k corx’esponding to 


and 


k'U = cj , for small k, 


a).." (k) = a cO ,for large k, 

u 


(29) 


(30j 


(30b 


with ^ more precise evaluation, based upon (28) instead of (29) , 

will give the same results (30) except with a slightly modified 
numerical value. 

C. Transformation of a Spectrum in 4),k-Spacc into a Spectrum in k-Spaci 
An intensity 


J i J w* 


(31) 


is obtained, by integrating (27) with respect to LO from-boto Co . If 
we perform tlie integrations with respect to dw ’ and d<J using (28) , vs/e 
reduce (31) to 



-11- 



j'j'dtj dk' y, A ( 


(32a) 

A - 



(32b) 


respectively. A comparison of the right hand sides finds the desired 
dispersion relation: 




Of, T/ny 

^ wv- 


C33) 


which serves to transform a function in c) ,k-space into a function in 
k-space, or reciprocally. 

This dispersion relation is no^^; used to transform (22b) into 



. V. SPECTML STRUCTURE OF PRESSURE FLUCTUATIONS 
A- Frozen Turbulence 

In frozen turbulence, the free streaming at velocity ”17 predominates, 

i.e. 


reducing (29) to 


]t0'- k'-U I » dO , 

' l— I JJ) 

~ ^ ^ ' 


. and (3R) to 
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(35) 


with 


^(m) - Jciij'JcL^ A(^'/k'c^ <r("cj'-AV^^ 
-/ 

= fJc(^A(^M) 


f ■ M + J~^r^ 

2A7 \ 2 /- M 

to x'epresent the effect of Mach nximber M = u/c in a free stream of 

velocity TT. The function (36) is found to he independent of k’ . 

Note that the available domain of k’ -integration corresponds to the 

limits (0,k) in the k’ -integration, because the integrations refer to the 

truncated functions (]<») and (k’)| ^ • Tlius , 

— *-v^ V*. 

by means of spherical eooi’dinates , we have v\/ritten 


(36) 


J cik^ ... - J dk' ^rr 


-I 

in (35), where U = k’*u/l<’u is the co-latitude, 

! •X<* WM 

In tlie inertia subrange (14b) , v^;e can calculate (15) , and by 
substitution, reduce (34) to 

f’dk' 


(37) 


or 


ii(i^ ) = jf/ 


(38) 
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in terms of j|0O from (16b). 

B. Fluid Pressure in Incompressible Turbulence, i.e. c =oo, or M=0 
With c = Co , we have ^ (tj Vk/c) = 1, and reduce (34) to 

after integration with respect to d(j* using (28) . Further calculations 
which need not be elaborated give the same result as (38) , with 


0~j^(M=0) = 1 


(‘ 10 ) 


C. o -Spectrum in the Absence of Free Stream u 

The dispersion time (29) ,with u = 0, which peaks at (30) ,can be used 
to calculate (39) - We obtain the following intensity 

>= („i) 

Co> 


and the spectrum 


where A is defined by (5) . Since the Mach nujiibcr 

^ 1^ ^ ^ 

may be considered as negligibly small, wc can reduce (M2) to 

cj-\ 


(M2) 


C'I3) 


The calculations which lead to (41) involve the use of (30) , a 
change of variables from k’ into according to (25) , and an 

integration with respect to . The details are omitted. The niunerical 


coefficients b, , b depend on (£ as ; 

t- 
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P M “3 

•3 


( 44 ) 


The power law O , as predicted by ('13) , can find agreement with 
experiments of Gorshkov CO . 

D. 0) -Spectrum in Non-frozen Turbulence with Free Stream u 
The spectrum, as determined by (3'l) , is governed by a product of tv\/o 
factors 


(4,’A’c) ^ (k») I and (k») , 


015) 


and is weighed by a dispersion time Tl . From (27) , (29) and (30) , 

this dispersion time is seen to be bi-modal: IVlien the first factor 
predominates wliile tlie second factor is insignificant at small k’ , 
has the property (30a) ; on the other hand, when the second factor predominal 
while the first factor falls off at large k’ , TL has the proper>ty (30b) . Thu 
the integration with respect to dJe’ in (34) , which is weighed by T in the 
manner described above, amounts to 

k''^A(o'A'c) Af%'jiy s(o'-k'-u ) 






017) 


The enstrophy function in (47) is 




= 


q, cj' 

Fi 


(48a) 


^ CJ 
u 


2 . 



-15- 


with 



and the remaining integral performed in spherical coordinates is 


= zA c ^{m) co' 


"A 


witli 


"-ti ^ ■ 

By collecting (47) and (M8a) and substituting into (47) and (34) , 
we find an intensity 


A </^'^> -- J ^ )j C(c^' C^' 

and a spectrum 


60 


- 5/3 


I 2 


h . i/i 


- 5/j 


The power law -5/3, as predicted by (52), agrees witli experiments bl 


(48b) 

(49) 

(50) 

(51) 

(52) 
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APPENDIX: EMISSION AND SCATTERING OF SOUND BY TURBULENCE 

The propagation of sound, as represented by the variables 
l> = density, p^=pressure and u^= velocity, can be obtained by an operator 
cC which selects longitudinal mode, as may be driven by a turbulent source. 
The tx’ansversal modes are then selected by a complementary operator 
1 equations of continuity and momentum of the sum of 

longitudinal and transversal modes are: 




(Al) 


(A 2) 


the effects of viscosity being neglected. 

Upon differentiating C^-l) and (A2) witli respect to t and x, and by 
subsequently taking the difference, we obtain 




^ y- f )■ ■ 

by applying the operators ^ and j[^ , we obtain: 


CA3) 


Vf^ 


-7^ A ^ 




y 




L T T Al 


(A'lb) 
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Note that the density in the transverse mode is simply r = constant. It 
is obvious that the sum of (A4a) and (A4b) recovers (A3) . 

We assume that turbulence is completely transversal on account of 
its incompressibility, and is not much modified by the propagation of 
longitudinal modes, and that the sound, vvhich is represented by 
longitudinal modes, lias a linear propagation. Then we reduce (A4a) to 




(ASa) 




(A5b) 


f 2 T T 

V. -u u . = r represents a turbulent source for the 

' I L 

emission of sound, and the terms linear in and u represent 

scattering. 

By neglecting scatterings , and linearizing tlie longitudinal modes, 
we obtain the Lighthill equation in the foi’in 
2- nL -1-1 L. ^ 




or 


(c^ V - P 


^ A, 


CAG) 


if the equation of state 


dp/d/3 = G 


f 


with the speed of sound c is used. 





Propagation of Light Through Random Refractive Index 
Fluctuations in Non-Frozen Turbulence 

C- M. Tchen 
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ABSTRACT 

A nonlinear dispersion relation for the space-time 
transformation in light propagation through a turbulent medium 
is developed. The derivation is based upon a nonlinear propa- 
gator, calculating the perturbed trajectory along which an 
optical field is transported by turbulent motions. The non- 
linear theory of dispersion relation is developed without the 
hypothesis of frozen turbulence. The frequency spectrum thus 
found agrees better with experiments than the one based on the 
hypothesis of frozen turbulence. The remote sensing by two 
laser beams can measure simultaneously the mean cross-wind and 
the eddy diffusivity in this method. The dispersion relation 
without the hypothesis of frozen turbulence finds applications 
to other remote sensing techniques, e.g. remote sensing by two 
receivers, and by multi-wavelengths with or without absorption 
by aerosols, and methods of reduction of turbulence effects in 


remote sensing. 



I . INTRODUCTION 

Recent progress of laser technology creates a great 
interest in optical turbulence, dealing with the coupling be- 
tween laser irradiation and turbulence in the atmosphere. We 
shall discuss tlie following: 

(a) laser light propagation in the turbulent atmosphere, 

(b) remote sensing of the turbulent characteristics and the 
structure of the atmosphere, 

(c) image retrieval, 

(d) propagation through aerosols. 

Since the wave field is coupled to the turbulent fluctuations 
of the refractive index, a parallel effort for investigating 
the statistics of wave field and atmospheric turbulence becomes 
necessary. 

In the statistics of light waves, the knowledge of tur- 
bulence and of the turbulent characteristics of the atmosphere 
has been meager and limited to the assumptions of validity of 
Kolmogoroff ' s spectrum and of Taylor's hypothesis,^ although 
the turbulent motions in the atmosphere are of large scales 
and will clearly deviate from these conditions. Even in small 
scale turbulence, the mean wind may be weak or there may exist 
other circumstances where these approximations are not valid. 
These limitations hamper a correct interpretation of optical 



3 


data for retrieving the atmospheric characteristics. 

The reasons why these assumptions have been made are 
evident. The assumption of a Kolmogoroff spectrum simplifies 
the statistics of wave propagation to the problem of a sto- 
chastic process with one single parameter. The assumption of 
Taylor's hypothesis or of frozen turbulence relieves us from 
the difficulty of a space-time transformation that is generally 
valid for frozen and non-frozen turbulence conditions. The 

importance of the general space-time transformation has been 
2 

pointed out, and there have been suggestions of substituting 

the standard deviation of wind velocity fluctuations for the 

3 

mean wind velocity, if the latter is too weak or absent. 

The linear theory of light propagation produces a light 
scintillation (P linearly proportional to the structure 

coefficient (? of the refractive index fluctuation 'K, , 

/fX, I 

as an index of the intensity of turbulence in the inertia sub- 
range. Experiments show that this linear relation ceases to 
2. 

exist when (T ^ O-J , a limit easily surpassed at a large 
optical path-length [_ in space probing. Under this circum- 
stance, a nonlinear theory of light propagation is necessary. 

On the one hand, a nonlinearity can appear when the light 
field E is sufficiently intense as to heat the turbulent 

} E j 


medium, so that 



4 


Then the equation of propagation of E -field will contain the 
nonlinear terra 

We 

whose role is to produce a self-focusing of laser beam and a 
therraal blooming. On the other hand, the propagation of a weak 
E -field in a refractive medium is governed by an equation of 
propagation linear in and E . But the presence of the 
coupling term 

requires a nonlinear statistical treatment for closing the 

statistical hierarchy entailing moments of ever increasing 

high orders, even when the statistical characteristics of 

are given. The literature treats this coupling under the name 

4 

of "saturation." Experimental devices using differential 
transmitter and receiver apertures can suppress the sensitivity 
of large-scale turbulence and thereby correct the effect of 

5 

saturation. Nevertheless, the optical effects from large 
eddies and the resulting saturation remain to be an important 
task of measurements for the purpose of understanding the 
large-scale structure and characteristics of. the atmosphere. 

In the following, we shall investigate the space-time 
transformation by considering the nonlinear perturbation of the 
trajectory along which a scalar quantity is transported by 
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turbulence. This nonlinear dispersion relation plays an im- 
portant role in the statistics of optical field and remote 
sensing. The new parameters which appear in the formula of 
transformation provides us with a means of remote sensing not 
only the mean wind, but also the turbulent characteristics. 
Illustrations are given for: occultation, multi-wavelength 

laser with and without absorption, cross-wind measurements, 
and frequency distributions. Comparisons with measurements 
are made, and show the importance of the nonlinear 
dispersion. Large eddies having a spectrum different from 
the Kolmogoroff spectrum are also considered. The space-time 
transformation which enters in the problem of propagations will 
also play an important role in "image retrieval." 
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2. PROPAGATION OF WEAK LASER IRRADIATION 
2.1. Wave Equation 

The propagation of light is governed by the following 
Maxwell equations: 


„ pi , 

U I 

V K ti - -L — 

•■V 


( 2 - 1 ) 


/V * ^ 

= o, 1> = £ e 


V- fy = 


/V 


where UU) is the electric field, H(t^) is the 

magnetic field, £, is the dielectric coefficient, and C is 
the speed of light. 

By eliminating f~j , noting that 


Vx7/ E --V^E -t V V-E 

_ _ v'-f ~v(b- E ) , 


(2-2a) 

(2-2b) 


and 


2 _ 

--t 




€ 


y ~ 

B 


T) 


t' 


(2-3) 
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and assuming that the atmospheric inhomogeneities vary slowly 
with time, i.e. 


= £(^) 


we obtain 


c 

<7^E 1 - v(e- viLi) ¥ ^ € I -o^ 

o 

c oSt 


or 


( 2 - 4 ) 


( 2 - 5 ) 


when we let 


E 


Here 


k 

r> - » ^ 


(2-6) 


is the optical wavenumber, ^ is the unperturbed dielectric 
coefficient, and the perturbation is 




( 2 - 7 ) 


It can be assumed that the scale of atmospheric turbulence 
is much larger than the optical wave length, so that the de- 
polarization as represented by E- in (2b) , or 

i*4 V4* 

in (2a), can be neglected, obtaining the wave equation for the 
scalar variable f in the form: 

+(C/Qt E ^ o . 


( 2 - 8 ) 
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For the dielectric coefficient, we can write 


giving 

k + (K 

^/ = ^ »*7 ( 

where H is the refractive index, with an unperturbed value 

4t and a perturbation . If the effect of yv, on the 

P 1 / 

propagation of is weak, and so is the heating of the at- 

mosphere by the laser, we can neglect the quadratic term, and 
write 

4 = £, 3 , ( 2 - 11 ) 

reducing (8) to 

B 4- ( 2 - 12 ) 

For the sake of abbreviation and without loss of generality, 
we shall write (12) as 

^ ^ f /t^ 2 E — O ^ (2-13) 


^ ^ ^ (2-9) 


> 


■tv. 


(2-10) 


- <0; - 


with ^ ! 

0 I 


and 


jx = ox, ! 
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2.2. Rytov's Method of Solution 

The problem of single scattering can be best explained 
from the wave equation (13) as a point of departure. The 
Rytov method introduces a function 




(2-14) 


with an unperturbed component and a perturbation ^ , 

'o / 


such that 


c/ I P/ I 


(2-15) 


This transforms (13) into the following linearized system: 


y'-A +zvf^vf i-Z k^\ = 0 

* Wl bv> 


(2-16a) 


(2-16b) 


Another change of variable 

-t 


i- 


V -t 


(2-17) 


transforms ( 16b) into the differential equation 

2-2 ^ 

V j- k -h z k E = O 

0 P / 


(2-18) 


which gives the solution: 
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(2-19) 


or 








G" (x') . 
7 ' ' o ' 


( 2 - 20 ) 


The integration extends to the whole available domain. 

For the sake of simplification of writing, we introduce 
a Green's function 

^ i \ . u k U 1 

e (2-21) 

^^1x1 

and a kernel 

{{x, x') = (^l X-X'l) E:/x')/e(^) , (2-22) 

simplifying (20) to the form: 


f.U) ~ 

' ^ tv V4^ ^ 


(2-23) 


If the scale of turbulence in the K direction is much 
larger than the scales in the -directions, the contribution 

to f'l in (23) comes mostly from small values of I ^2 " | • 

w 

This permits one to approximate | X j by 






(2-24) 



or to approximate (21) by 


<rOrri)“ 




lij. 


x~x -i- 


J. 


X -X' 


and (22) by 

transforming the solution (23) into 



We can write 


'f - A i- 


(2-25) 


(2-26) 


(2-27) 


(2-28) 


as decomposed into an unperturbed value 

t = A ^ 


(2-29) 


and a perturbation 


from 


t, = ^ ^ ^ 






(2-30) 


with 
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Here is called log-amplitude fluctuation, and Sy is 

called phase fluctuation. 

Upon identifying (27) with (30), we find 


S,G) 




(2-31a) 


(2-31b) 


where 


'ty - 4^ -h L 4. 


(2-32) 


2.3. Solution for Plane Wave 


We consider a plane incident wave 

tkK 


^ok) = 


(2-33) 


transforming (27) and (26) into 

• X 


0 VM 


(2-34) 






I ,'ck 


2-tc ;(-x' ' zC%-)^') \ 


(2-35) 
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1 


J 

J 

J 

;j 


J 

] 


The convolution form, as present in (34) , suggests a 
Fourier transform in the transverse direction. Thus by writing 


'VI 




(^y i J / 


with 


5^ ^ (j'j) - A " ) 


we transform (34) into 


with 


(2-36) 


(2-37) 


{^x'j =: ik^ 1^- V 

Since /l, is complex, by (32), we find 

j k] k ) 

k) KJk 


(2-38) 


(2-39) 


(2-40) 


( 1 



with 
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3. STATISTICS OF FIELD FLUCTUATIONS: CORRELATION AND 

SCINTILLATION 
3.1. Plane Wave 

The correlation function of "X -fluctuations is defined 
by the mean product 

+ f fy ( 3 - 1 ) 

of the X -fluctuations at the two points 



w»-» 


on the plane transversal to the direction K of propagation. 


Here and in the following we omit the subscript " "in X 


and . 


By the Fourier formula, we can write 

V?) -/•'i l(i) . 

with the aid of a spectral function 


■/ 


(3-2) 


■ < 3-3 

Here J is the infinitesimal surface element, 

and ?.)( is the interval of Fourier truncation. 

For the calculation of (3) , we write 

~J (^r) ^ ^ 


(3-4) 



at K- L 


from (2-40) . We can formulate 




^y<(ijk)^(Ly-k)y =- J [zn)‘^yLj.L-x'kjlJ^L~ J 


X ~yy 

I 


(3-5) 


For large /_ as compared to the scale of turbulence in 

the X -direction, we can replace the integrations with 
respect to X' X '' by integrations with respect to 


Xp , K, X'- 


as 




0 -Oo 


The integration with respect to Xy takes the limits 
^-Cto^OcJ approximately. The kernel varies slowly 

with X' and x" , so that we can approximate 

-X' k) = 4^ (L-x'j k ) . 


Hence we can reduce (5) to 


'L 


- J NAji-K'ijJ 


(3-6) 
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The last integral is 
-^^o 


' ~fiO ' 


(3-8) 


and can be expressed in terms of the spectral function 
of 'Tu -fluctuations, so that (3-7) becomes 


L 




(3-9) 


It follows, from (2) : 




(3-10) 


The integration with respect to ^ can be written in 
polar coordinate, as 


oJl, f/ ~ 4y hikk 2,-ri 

A Vf. J V 

0 ''o 0 



''In I ^ 


= y it 'j 4k k z ^ V] 


= j 4 t'Jikk Z(k/’)^(/‘J '^^(\k) , 


(3-lla) 


'‘o 


or as 
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00 

k l(kf) i(k) %ik], 

with the use of a filter function 



W y y> / 


Here 


7 =• 





We note that ^ 

represents the Fresnel zone size, with A=2’7r/k^- 


By definition, the structure function is 

^ (^/() ‘^f) -/ ( )] y 

K4 

= - K,(uj] ■ 


(3-llb) 


(3-12a) 

(3-12b ) 
(3-12 c) 


(3-13) 


(3-14) 


(3-15) 
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Since is given by (11b) , and ^ is given by 




= ZTt^k^L J dk k $Jk) f^(k) , 


(3- 


we can rewrite (15) as 


f u k $jk) 

•V, Jq 

The statistical characteristics of S-fluctuations are 


obtained by considering 


i. 


instead of 




‘A 


Analogous 


calculations yield the following formulas: 


k fjkl 

0 


"O 

stO 


(3- 


(3- 




with a filter function 


l(k) = I + -^^ 7 ^ 

r 

- [ 7^ / 

- I ^ fv^ 7 > / . 


(3- 


16) 

17) 

•18) 

19) 

20) 

•21a) 


(3-21b) 



The most usual spectrum is the Kolmogoroff spectrum 


ijk) = 0.0^3 C 


(3-22) 


It is valid in the inertia subrange 

< k < . 

Because of the presence of the filter function -fjk) ■ 


(3-23) 


the 


contribution of the product 




to the integral 


will depend on the value of relative to 


and k 

V 

. In this connection, we 

shall 

consider the 

following 

three regions: 




(a) 

geometric optical region 

K 

<' 

(3-24a) 

(b) 

diffraction region /C 

O 


< 

(3-24b) 

(c) 

large scale turbulence region 

kp < K • 
r D 

(3-24c) 

The three 

regions are illustrated in 

Fig. 

la-c . 



In the region (a) , the tail of the Kolmogoroff spectrum 
near the cutoff at 

contributes most significantly to . By (12b) , we 


have 


- ‘'A, 


2.'7U^k^ L J (Ak k k ) 


0 
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with 


i. ^ ar" 

7 3 ; 

In the region (b) , the Kolmogoroff law contributes 

(■)(^) = O.C 5 J vw'^k^LjcLk k ~ , 

O Z" 


(3-26) 


by (16) and (22) . Since the integral is 

-573-r“, -8/j 


0 

we find 


y^;7~ 




(3-27) 


(3-28) 


In the region (c) , the contribution of the Kolmogoroff 
law (22), that is restricted to the inertia subrange (23), has 
the form 

A 


f 4k k $Jk) ^ 
VC 


(3-29) 


or 


k 


Lrk'L f Ak k $JkJ 




(3-30) 
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with 



The variance (30) has a milder dependence on L than (25) 
and (28) . 

Often the turbulent medium has a spectrum not merely- 
restricted to the inertia subrange (23) , but possesses in 
addition a spectrum 

in the subrange 

k ^ K, 

c£ larger-scale fluctuations concomitant with (12b) . This 
spectrum will contribute a variance 


-31) 


- 32 ) 


with 
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The parameter C depends on C and Yk - The spectrum 

•h. V- 

(31) takes its origin from an inhomogeneous 


3.2. Spherical Wave 

The spherical wave 




) 


I 


‘t-Trjxl 

✓A 

' ' - f «.( 




X + 


(3-33) 


f 7T X ‘ ' ' 2L X 

has a form similar to the Green function Q~(>^ ) (2-21) 

and (2-25) . 

Upon substituting (33) and (2-25) into (2-26) , we get 


'i'J = 




— -cxi i y=X'/x <3-34) 

By comparing (2-34) and (34), we note that (2-34) for the 
plane wave can be converted into (34) for the spherical wave 
by multiplying all variables except in (2-34) by y* . 

In this manner, (2-41) is converted into 




S (2u:) ^ ) 

to^ 1 -kCl-x ') = V 


(3-35) 
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and (3-lla) is converted into 


k jjnf)ia)^'%;k) . ( 3 - 36 ) 

^0 0 ^ 

By letting ^ (36), we obtain the variance 


(3-37) 


When the Kolmogoroff spectrum is applied, we find 

(Tj ^ = O - OlS (2x)^kJ'J dy, ' C^(r) J dk k Y ^ K(L X') 


(3-38) 


The last integral is 


r I r(L -*'il ~ ^ 

hy ^^7 ' 


so that (38) becomes 




(L-xO 


L L 


t/i 


, b( =. 0-S(}3 


(3-39) 


(3-40) 


/«-s ^ 

If = constant, (40) reduces to 

•71 

/ , _ 2 / l/C, I n/^ 

^ ^3 = 0-/2^. ^ 


(3-41) 
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having a form similar to the plane wave case (28) with a dif- 
ferent numerical coefficient. 

The formulas (40) and (41) are valid for the diffraction 


region (24b) . 
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4. REMOTE SENSING OF A PLANETARY ATMOSPHERE BY THE 
METHOD OF SCINTILLATION 

Consider a planet of radius ^ . The light propagates 
from a transmitter ~J~ (a space satellite) , placed at X = O, 
and is measured by a detector placed at x''— [_ , See 

Fig. 2. The point Q on the line of sight is situated at 

and is at the shortest distance ^ from the center 
P of the planet. A variable point M on the line of sight 
is situated at x " . 

We can assume 


and write 


L « ^ ^ 


A = 


+ (x'-L,/ 




(4-1) 





(4-2) 


The remote sensing by the light propagation from a 
space probe can measure the structure of the atmospheric en- 
vironment. Applications can be found in weather prediction, 
climate modification, pollution studies, storm and clear air 


turbulence warnings. 



1 
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We shall discuss the following case of spherical wave 


propagation. 


If (2 depends on the height /L- from the surface 

'IV 


of the planet, as 


^2 Z - H 


(4-3) 


where j-j is constant, we can substitute (2) and obtain 
Z. I * ki^(^'-L,) ] - 


(4-4) 


A further substitution of (4) into (3-40) finds 


7/^ ^2 


^r}-,k'cU 


4 “ l”r“l ‘•t [- I . 


(4-5) 


In view of large ^ , we can approximate the integral 

H 




L, 


= iz; 

ku \ 


(4-6) 
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and obtain 






/ / 1 

(L,L^/l) ^ a^0-Si3. 


( 4 - 7 ) 


This result is valid in the diffraction region. We find a 

Lj and . 


reciprocity in 
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5. STATISTICS OF FIELD FLUCTUATIONS: FREQUENCY 

SPECTRUM IN A CROSS-WIND 
5.1, Frequency Spectrum 

We consider a plane wave propagating along the 
X -direction. See Fig. 3. Two receivers are located at 
2^ and ^ cross-wind of velocity U, - 

vv% kv* 

and measure 

dn \+?) , (5-1) 

to form' the correlation function 




a 






(5-2) 


with fixed ^ and variable 

We consider a wind of velocity 

H — U + M, ^ 

.V. 

consisting of a constant average velocity M, and a fluctu- 


ation q 

Since the pattern J st time 'ttT- 

position must have come from time "t and position 


and 


- UZ -JUxJ 


(5-3) 
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we can write 




(5-4) 


in 


on 


the notation of a propagator , or evoluti 

operator, w>iich indicates that the transport of ^ from it" 
to 't-f'C' is made along the trajectory perturbed by turbu- 
lence- Thus (2) becomes 




= /[f; - 







(5-5) 


The derivation follows that of (3-llb) , with an additional 
factor ^ representing the effect of Z(x) in the 

perturbed trajectory. The detervaination of will be 

given in Sec. 7. Note that the case with 


^(k^Tj = I y C>A A(~^) - 0 


(5-5) 


is equivalent to frozen turbulence. 
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By letting ^ 0 in (5) we find the 
temporal correlation, as follows: 


B^(t, 




(5-7) 


In wave propagation, a frequency spectrum is defined 


as 


00 


-bo 


(5-8) 


and becomes 


W^(wj ^ ^lo^k^Lf 4k k jjku-^^(k-) , 

t5> 


or 




hlYoj =9%\^L fdk k 


I 


(5-9a) 


(5-9b) 


Note that this definition differs from the usual one 


n / _ tCOT" 

-bo 
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with 


Mh -j h e ]J_k u-cj ^(k,z) 


(5-10) 


In particular, if the hypothesis (6) of frozen turbulence is 
valid, we can simplify (10) to 


- j (fz e ^ J (kuz) 


o, ka<t^. 


(5-11) 


and reduce (9b) to the well-known formula of frequency spectrum 
in frozen turbulence: 


u/u 


(5-12) 


5.2. Remote Sensing of Cross-Wind 

From phase fluctuation S, we can obtain a correlation 

function ^ ^ formula analogous to (5), if ^ 

5 ' ^ / I p( 

is replaced by / . The structure function 

15 




(5-13) 


i 
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takes the form analogous to (3-20) . 
The time derivative is 






1=0 


= j ‘tk k 




(5-14a) 


In the hypothesis of frozen turbulence 
degenerates to the well-known formula 


i'?=0 ' 


= - Ik k\(k)Mj J(if) ( 5 -i'ib) 

T=o '>0 ' 

for the determination of the mean cross wind U . On the 
other hand, in the present formulation where this hypothesis 
is not necessary, the general formula (14) can determine both 
the mean wind U. and the turbulent characteristics as repre- 
sented by . As one important characteristic property, 

the eddy viscosity can be measured directly, without going 
through the intermediary of eddy stress and velocity gradient. 
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6. REMOTE SENSING BY FREQUENCY SPECTRUM FROM TWO WAVES 


] 


OF DIFFERENT OPERATING FREQUENCIES 
6,1. Without Absorption 

We consider now the propagation of two waves of different 
operating frequencies 


Kl ' ' 


( 6 - 1 ) 


] 


with £^=1 / and measure the log-amplitude fluctuations 


and their correlation function: 


a(kzX,i)) 


L' > T = t-% . 


(6-3) 


According to (5-8) , the frequency spectrum is 




(6-4) 


- 6o 


Calculations find 




>2 


'0 ^0 


X/ 


(6-5) 


1 



35 








where 


cr/z \/ X2 ' 


^ i> ^ ^ 

*^o ‘‘o I 




oZ 


are defined by (2-41)^ and 


Y= x/L. 

- I, 


for a spherical wave 
for a plane wave. 


If we write 


(6-6) 


(6-7) 


H( - Q - 6 Qv 


(6-8) 


into two parts and , called cospectrum and quadrature 

spectrum, respectively, a coherence can be defined by 


Q^Xirkn'^) 


(6-9) 


This coherence has been measured, see Fig. 3. 



) 
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] 


6.2. With Absorption 

7 

Kjelaas et al. made an experimental investigation on 
laser absorption at two wave lengths, and showed that per- 
formance degradation due to scintillation can be reduced by 
exploiting the coherence between two fluctuations. The 
analysis is based on the propagation of two waves, described 
in Subsec. 6.1, but extended to include the absorption. 

The laser emits radiation alternatively at two optical 

with absorption at /?. 

only, giving the intensities at the receiver 
. ^ “ ^ 0 / ^ \ I 

-OrL 


wave numbers and with absorption at /? , 

0/ ^2. ^ o/ 




cj: 'k 


02 


(6-9a) 

(6-9b) 


in a laminar medium, where is the absorption coefficient 

of the pollutant gas of concentration C , while OJ are 

losses due to aerosols and other gases, and can be assumed of 
having a negligible difference 


(T, -(T^ ^ 0 


It follows from (9a) and (9b) 


C = 


o^L 


- Ih 


Zo 


I 


/o 


Z 

oLL 




(6-10) 
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or 

since has no absorption. The expressions (10) and (11) 

are valid for perturbed concentration as well as for perturba- 
tion. In the following, we shall consider a, 
perturbations or turbulent fluctuations. 

We can form the correlation function of (11) and obtain 

\(v] = (c(i) c(i+T^)y 

= kAQ' ([^^(t^^J-7/V] 7M]) 

(6-12) 

where 

is the correlation function of J)( , with 

A frequency spectrum in the form 

-t'lOT 


(6-13) 



(6-14) 
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or 


- ^ f ‘t-c ‘'^t^aJ, 

'^-OO 


(6-15) 


can be written in analogy to (5-8) . By means of these defini- 
tions, we can transform (12) into 




h/, f iV L^) -2. Ill ^ uT 

Xz2 /!! 


(6-16) 


The result is shown in Fig. 4. The spectrum decreases as 
Q and ^z/ X ! increases to approach unity. 

It remains the problem of derivation of the spectral 
structure of c-f luctuations from the spectral structure of 
refractive index fluctuations. The problem involves again the 
transformation between the frequency and the wave number spaces 
in a non-frozen turbulence. 
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7. SPACE-TIME TRANSFORMATION IN NON- FROZEN TURBULENCE 
7.1. Transformation of Eulerian Correlations in Time and Space 
The analysis of frequency spectrum has usually been made 
using the hypothesis of frozen turbulence, in which the pattern 
of refractive index is assumed frozen during its transport by 
a constant velocity LL . In the following we shall analyze 
the frequency spectrum without this hypothesis. To this end, 
define an Eulerian spatial correlation 



and an Eulerian temporal correlation 


5 (a(t, 

If the dynamics of ^ ^j_ J is governed by 

- 0 , 

with a differential operator 

L - 


(7-lb) 


(7-2) 


(7-3) 


of transport by a wind velocity U ^ U U. y consisting 

of a constant mean velocity U and a fluctuation u , we 

vV' 

can conceive that the turbulent pattern, as described by 


'/(tf-C, ) 


(7-4) 
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and occupying the position at time ~t -fZ' , must have 

originated from the position time ^ in a 

pattern described by 


?([t, fj, -Mr)] ^ 'U(t+X,t) 

The propagator is an evolution operator that 

represents the transition from iM to 't-t'C- along a tra- 
jectory ^ perturbed by the fluctuating velocity from (2) 
and (3) . The trajectory Jc obeys the dynamical equation 


. 0^1 4Mt 


✓1 

lA. 


(7-6) 


By equating (4) and (5) and substituting into (2), we 


have 





In Fourier form, we can write 




(7-7) 


(7-8a) 


and 


V(ttMj Mt) = ^ [t, -)Urj] 


-If I !L -ik-ihj 

-JjdudkA, ~~ ^ 
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and find the correlation (7) in the form: 



(7-9) 


Here the conditions of statistical stationarity and homogeneity 
have been used, and the hypothesis of independence 



by (3-3). 


By following a procedure analogous to that which has 
led to (3-lla) , we transform (9) into 
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By identifying with (5-5) and noting that 

f^(k) 

we determine 



(7-12) 


(7-13) 


It is thus easy to verify (5-6) in a frozen-turbulence. Thus 
the distinction between the perturbed trajectory in non-frozen 
turbulence and the free trajectory in frozen turbulence lies 
in 


^ j S - // respectively. 

We conclude that the Eulerian temporal correlation is thus 

obtained from the Eulerian spatial correlation that is 
represented by its spectral function . 


7.2. Calculation of the Perturbation Kernel 


The perturbation kernel 






■to 


(7-14) 



/ 
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can be calculated by means of a probability of transition 


yt j that governs the random path 


in a time IT / 


and is determined by the Fokker-Planck equation 


It 




(7-15) 


in isotropic turbulence. The eddy diffusivity K , which is 

'V 

characteristic of IL -turbulence, in accordance with the 
trajectory perturbation (40), takes the Kolmogorof f-Richardson 


form 




^ I 


(7-16) 


when the spectrum of velocity fluctuations obeys the Kolmogoroff 
law, with a rate of energy dissipation £ 


In Fourier space, (15) is 


where 




:,k) 


-IkJ 


i 


(7-17) 


(7-18) 


is the Fourier component. Noting that 




(7-19) 
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by the condition of normalization 

j<li [y(r,A ) = / 


we find the following solution 



and the perturbation kernel 





from (18), (19) and (21). 


(7-20) 


(7-21) 


(7-22) 


7.3. Frequency Spectrum 

Now we substitute (22) into (11) , to obtain 







(7-24) 


from (11) and (12) . 


The approximation assumes uniform 



and tc in the space. 
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By the Fourier transform in the definition (5-8) , we 
can transform the time correlation function (24) into the 
following frequency spectrum 


^ j ^ ^ ft ) fy (k) !<; K, u) 


(7-25) 


In an analogous manner, we derive the frequency spectrum of 
phase fluctuations 


5 

Here 





IS 


^4 U 


(7-26) 






-Ki' 




(7-27) 


/ 


0 


is a weight function. It takes the following two asymptotic 
values: 

(i) frozen turbulence with til » /C/g or k »k 






;? 6j 


o 


k i\. 


(7-28a) 

(7-28b) 


(ii) non-frozen turbulence, with ia. ■\< Kk 


k , or 


k<<k 


K4 V ^ ^ 




(7-29) 


The critical wave number ^ which separates the two 


cases 
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(i) and (ii) is 

The numerical coefficient 


(7-30) 


=0.9G (7-31) 

8—10 

has been determined by Tchen, using a kinetic method. 

Note that the integrands 

i 

and 


it J/V M(u,k, 



in (25) and (26) are functions of the following two dimension- 
less variables 



and 



(7-32) 


which characterize the three regions (3-24a) - (3-24c) of 
diffraction and distinguish between the frozen turbulence (i) 
and the non-frozen turbulence (ii) . 

For a spectrum lying in the diffraction region 

(3-24b) , we have 


(7-33) 
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The remaining parameter 




from (32) characterizes 


the following three cases relative to the degree of frozen 
turbulence: 

ijk^ - \ ^ << I ^ ^ » /. 

In the asymptotic limit of high frequency oO oo 


have: 


(A) 4 - 


A/ 6o 
Aj LO 


~ 2 . 


- 8/3 


corresponding to (ii) 
corresponding to (i), 


(7-34a) 

(7-34b) 


(B) ^ L 

u, r 


-2 


corresponding to (i). 


(7-35) 


corresponding to (ii) 


(7-36) 


^ y Vv^ 

In the asymptotic limit of low frequency ( O ), we 


have 


^Tt%'L J4U. ijk) f/k)h4.[u^o,hk;c:J, 


(7-37) 


t (u9 o] = J oik k 1/ 



IS 


/ ' ^ 


O 
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with 



from (25) , (26) and (27) . 


(7-38) 


i 


From the properties (3-12) and (3-21), it is seen that 
weighs toward small scales and weighs toward large 


scales, so that the frequency spectrum Ia/I is dominated by the 
frozen turbulence, while the frequency spectrum is dominated by the 
non-frozen turbulence. These physical features and the spectra 
(34a, b) have been observed in atmospheric environments. See 

TT- c 11 

Fig. 5. 


The time-space transformation in non-frozen turbulence 


is developed above by using the approximation of independence 

lO 

(10) . A correction has been estimated, and is found small 


in isotropic turbulence. 
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8. METHODS OF REDUCTION OF TURBULENCE EFFECTS IN 
REMOTE SENSING 


8.1. Principle of the Methods of Reduction of Turbulence Effects 
In remote sensing, an unperturbed field from a 

source or object at the transmitter plane ' is perturbed 
by turbulence, so that the field iyl(^ J at the receiver takes 
the form 




W = / Uf, r) ^ 

- 


( 8 - 1 ) 


where h ( ^ is a kernel function representative of the 

perturbation. For interpretation of measured data, the principle 
is : 


A. to solve for tC, , or some statistical characteristics 

of if kernel h can be found, 

B. to minimize the perturbation by means of measurements 
at two wavelengths. 

We shall describe below four methods under principle A, 
and one method under principle B. 


8.2. Speckle Interferometry 

This method was developed by Gezari, Labeyrie and 
12 

Stachnik. 

The perturbed field (1) is transformed into irradiance. 
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yielding a relation between the perturbed and unperturbed ir- 


radiances I, I , in the form of convolution integral 
o 


1(f) = Ijfj * ijif) / 

which is solved by Fourier transformation, yielding 


( 8 - 2 ) 









(8-3a) 


where H(j‘' ) and H(K) are called "spread function" and "transfer 
function," respectively, and K is wave number. 

Since a star can be considered as a point source, i.e. 
a -function, one can determine H(K) by pointing the receiver 

to a star, and calculate by 

from (3a) . 


13 


8.3. Intensity Correlation 

This method was developed by Korff. 

Instead of developing an equation for irradiances at 
one point, as with (2), it was suggested to consider the 
correlation function of irradiances at two points: 


CM- . 


(8-4) 
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A convolution is developed for correlations C^^(r) and 
C(r) at the source and receiver, in the form: 

c(k) = I C ■ (8-5) 

^ '] 

If is a slowly varying function, and H(K) is a 

rapidly varying function, (4) can be approximated by the 
relation 

C(KI wfdK'HiK'), 

li>, ( 8 - 6 ) 

which determines the spectral structure of C^(K). Then a 
subsequent inversion will yield the correlation at the 

source. 

8.4. Phase Retrieval 

The methods described in 8.2 and 8.3 do not provide the 

14 

phase information. Therefore Knox suggested the retrieval of the 
phase information by the following formula: 

- JJ^K) lXK^ iB-v 

8.5. Image Compensation 

15 

This method was developed by Deitz . Instead of one 
single processing of the field at the source through 

turbulence with kernel h to measure Uj at the receiver. 
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following the formula 


.00 







f-f2 


(8-8a) 


-iV 


as in (1), this image is subjected to a second processing through 

J 

a lens, or "imagery," with a known kernel Aj » to yield a 

second image at a focal distance, according to a similar 

r 

formula : 


^.(f) 


(8-8b) 


- Oc? 


The Fourier transform of both formulas (8a) and (8b), i.e. 


-- u^(k)1(kJ . 


(8-9a) 

(8-9b) 


reduces to 


lAo(kj-ijK) i(K) , 


( 8 - 10 ) 


An optimum reception is then obtained by a possible adjustment 


of A- such that 


ijK) 


rj 



-/ 


( 8 - 11 ) 


Similar processings can be made for irradiances in the 
place of fields, giving the relation 



This relation can be approximated by a simpler one 




(8-13) 


if 




varies slowly. 


8.6. Multi-Wavelength Coherence for Reduction of Turbulence 

Effects 

The method was developed by Kjelaas, Nordal and Bjerke- 
7 

strand, and has been discussed in Sec. 6. We give its principl 
in the following lines. 

In laser absorption at two wavelengths, the method by 
Kjelaas can show that performance degradation by turbulent 
scintillation is reduced by exploiting the coherence between 
the fluctuations at two wavelengths. 

The experiment finds fluctuations of concentration 
by measuring the log amplitude fluctuations 
wavelengths » separated by a time interval A , 

i .e. 




(8-14a) 


where is a constant, dependent on scintillation path and 

absorption coefficient. Similarly, 
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C- 


(t-f . 


(8-14b) 


From (14a) and (14b), the correlation function can be 
obtained, and is transformed into a spectral function for c', 
which takes the form 




where /^jj 

/(/ ^ / 


and 


W C(^) 


(8-15) 

are spectral functions of 


pc ^ ^ -fluctuations, and 


= C,/^i 


( 

'OIZ 


(8-16) 


is a co-spectrum dependent on A and * This method 

shows a reduction of performance degradation by minimizing 
\J\J [i^) > that is by maximizing * This is achieved 

by letting A and I ' method has been applied 

to atmospheric experiments. 



55 


9. EFFECTS OF ABSORPTION BY WATER VAPOR ON LIDAR 


16 


PERFORMANCE AT 10-6 jA. rr^ 

Recently experiments by Schwiesow and Calfee"^"^ have 
shown that the fluctuations of water vapor concentration are 
equally important, if not more important than the fluctuations 
of temperature in contributing to the refractive index fluctu- 
ations at • Experiments have shown that simple argu- 

ments of these new effects by water vapor are not sufficient, 
and would require a careful re-examination of the problem by 
an analysis of the propagation of light through the atmospheric 
turbulence. 

Most theories of propagation deal with a real refractive 
index fluctuation. But because the absorption by water vapor 
on lidar performance at 10.6 is important, we include a 
complex refractive index by writing 


'yi-h i'tny' 


(9-1) 


While a detailed analysis for a weak scintillation theory 
can be made in a separate task, we can predict the following 
log-amplitude fluctuations: 




(9-2) 
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IS 


where ^ is log-amplitude fluctuation, ^ 
optical wavenumber, 4 is optical wavelength, /L is path 
length, is the spectrum of atmospheric inhomogeneities. 

The index 


i = / J 


(9-3) 


denotes the contributions from 7iiv , and n couplings. 


The filter functions 




,) . are functions of h/k '< where 


% - (.K/l) 


(9-4) 


is the Fresnel wavenumber, and determine at what scales the 
spectrum should be effective. 

By writing 




(9-5) 


with 


fdi k iZkj/ aj , 


(9-6) 


we expect to find the following: 

(a) For i=l in nn- coupling, the filter function 
(kj determines that those scales in the Kolmogoroff inertia 
subrange are most effective, so that depends on 










2 Z 




t 


(9-7a) 
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and 



(9- 7b) 


in agreement with the weak scintillation theory without ab- 
sorption. Here ~J~ is temperature and ^ is density of water 
vapor. 


(b) For i=2 in coupling, the filter function 

determines that scales larger than the Kolmogoroff 
scales may not be negligible as in (a) . The results (7a) and 
(7b) are not valid, and will show additional contributions 


from 



(c) For i=3 in >n. 'Vt, coupling, we expect a law of 

which will consist of a mixture of (a) and (b) , by 

contributions from 

r n Jn 

-ff rf ' Tf- 

The large scale fluctuations of the refractive index, 
temperature and water vapor are not well studied in the litera- 


ture. 
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The complex nature (1) due to the presence of water 


vapor will find a pronounced effect in the reduction of 


turbulence effect when two optical wave lengths are used, as 

7 

suggested by Kjelaas. 
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10. DISCUSSIONS 

By means of a nonlinear analysis of the perturbed 
trajectory along which an optical variable is transported by 
turbulence, we have derived a space-time relation which greatly 
improves the statistics of wave propagation. As an illustra- 
tion, we have analyzed the frequency spectra for ^ and 
S-f luctuations, and generalized the formulas of remote sensings 
by extending to non-frozen turbulence (cross-wind determination, 
multiwavelength sensing) . 
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